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(57) Abstract: The invention provides a method for computer-implemented controlling of a processing apparatus having at least 
one independently controlled manipulated variable and at least one controlled variable responsive to the manupulated variable. The 
control approach uses a robust multi -variable controller which defines an expected variation in magnitude for each controlled vari- 
able as a respective function of each manipulated variable via use of a set of at least two models. The model set has a dynamic 
response inertial characteristic. The two models are derived from a Reference Model (the traditional model defined in Dynamic 
Matrix Control). The multi-model set further enables adaptation of the models and gains during real-time use. 



WO 02/03150 PCT/US01/18111 
MULTI-VARIABLE MATRIX PROCESS CONTROL 

This invention relates to a computer-implemented method and computer-based apparatus 
for adaptive control via a matrix of feed-forward relationships respective to a process or 
5 series of processes implemented in an apparatus such as a chemical manufacturing plant. 

Optimized operation of manufacturing systems is valued for providing benefits in 
profitability, productivity, environmental impact, and high product quality. The increasing 
capability of low cost computers to deliver resolution of complex control approaches has 

10 advanced optimized operation to incorporate techniques which could not have been 

economically explored even a few years ago. One of the techniques enabled at a relatively 
early stage of computer use in control was Dynamic Matrix Control, a form of feed- 
forward control based upon a method where outputs, or controlled variables, are predicted 
to move in the context of known control settings and current data. Feed-forward 

15 methodologies are, in many cases, superior to feedback methodologies which wait until 
process disturbances have actually changed the controlled variables before controller 
action is taken. Indeed, an ideal controller provides both feed-forward and feedback 
action in sufficient capability to achieve optimal operation. 

20 Dynamic Matrix Control is discussed in US Patent 4,349,869 for a DYNAMIC MATRIX 
CONTROL METHOD which issued on September 14, 1982 to David M. Prett, Brian L. 
Ramaker, and Charles R. Cutler. Dynamic Matrix Control has helped in solving control 
issues related to limitations in future controller response in the context of a decision at a 
particular time, accommodation of the full set of conditions in a system being controlled, 

25 complexity in multiple influences, and non-linear impacts respective to disturbance. 

But, in its classic form, Dynamic Matrix Control (DMC) also has presented challenges in 
dealing with process stability during transients or at different turndowns or deviations 
from the fluid flow throughput used in model definition and tuning. The need for stability 
30 in use of DMC in all operational contexts has been a driving force in the development of 
the present invention, which provides beneficial features in stability, robustness, and 
adaptability to substantially enhance the traditional DMC approach. 
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The invention provides a method for computer-implemented controlling of a processing 
apparatus having at least one independently controlled manipulated variable and at least 
one controlled variable responsive to the manipulated variable, characterized by the steps 
of: 

5 

defining an expected variation in magnitude for each controlled variable as a respective 
function of each manipulated variable from a set of at least two models, the model set 
having a dynamic response inertial characteristic; and 

10 implementing a change defined from the model set to modify the manipulated variable in 
the processing apparatus. 

In further detail, the invention provides a computer-implemented system for controlling 
the operation of a processing apparatus having at least one independently controlled 
15 manipulated variable and at least one controlled variable responsive to said manipulated 
variable, characterized by: 

a set of at least two models for defining an expected variation in magnitude for each 
controlled variable as a respective function of each manipulated variable, said model set 
20 having a dynamic response inertial characteristic; and 

means for implementing a change defined from said model set to modify said manipulated 
variable in said processing apparatus. 

25 The invention also provides a computer-implemented system for controlling the operation 
of a processing apparatus having at least one independently controlled manipulated 
variable and at least one controlled variable responsive to at least one said manipulated 
variable and further responsive to process disturbances originating independently of said 
manipulated variable, characterized by: 



30 



means for measuring the magnitude of each controlled variable; 
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a set of at least two models for defining an expected variation in magnitude for each 
controlled variable as a respective function of a manipulated variable disturbance instance 
in each manipulated variable; 

5 means for determining, from said set of models, an estimated modeling error value in 
interaction of one said controlled variable with all said manipulated variables in achieving 
a manipulated variable modification; 

means for determining an estimated process disturbance value from said controlled 
10 variable magnitude and said estimated modeling error value; 

means for defining a first portion of desired change in the present value of at least one 
manipulated variable from said estimated modeling error value; 

15 means for defining a second portion of desired change in the present value of at least one 
manipulated variable from said estimated process disturbance value; and 

means for implementing said first and second desired change portions to modify said 
manipulated variable. 

20 

The invention further provides a computer-implemented system for controlling the 
operation of an apparatus having at least one independently controlled manipulated 
variable and at least one controlled variable responsive to at least one said independently 
controlled manipulated variable, characterized by: 

25 

means for defining a set of consecutive discrete intervals of time in a time-dependent 
function; 



means for introducing, in each independently controlled manipulated variable, a 
30 manipulated variable disturbance instance of predefined magnitude, said disturbance 
instance prompting a response in each controlled variable; 



means for measuring the magnitude of each controlled variable; 
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a controller, 

means for defining, respective to each response, at least one time-dependent functional 
5 characterization of said controlled variable magnitude over said set of consecutive discrete 
intervals of time on a time-axis, said functional characterization having a zero-time time- 
axis attribute, a maximum-time time-axis attribute, a dead-time time-axis attribute, a 
response gain attribute, a ramp-rate attribute, a steady-state time-axis attribute, a 
curvilinear portion disposed between the time-axis position of said dead-time time-axis 
10 attribute and the time-axis position of said steady-state time-axis attribute, a ramped 

portion disposed between the time-axis position of said dead-time time-axis attribute and 
the time-axis position of said maximum-time time-axis attribute, and a homaloidal portion 
disposed between the time-axis position of said steady-state attribute and the time-axis 
position of said maximum-time time-axis attribute, each discrete time interval for one said 
15 response having the same time duration, said homaloidal portion having a zero value for 
an integrating controlled variable response, said ramped portion having a zero value for a 
non-integrating controlled variable response, and each functional characterization for one 
said response having identically valued zero-time time-axis attributes, response gain 
attributes, ramp-rate attributes, and maximum-time time-axis attributes; 

20 

means for acquiring a first said time-dependent functional characterization respective to 
the response from measuring an effected change in said magnitude of a controlled variable 
after introducing said disturbance instance, said first functional characterization having a 
first said dead-time time-axis attribute, a first said steady-state time-axis attribute, a first 
25 said curvilinear portion, a first said homaloidal portion, and a first said ramped portion 
having its functional derivative equivalent to said ramp-rate attribute at said maximum- 
time time-axis attribute; 

means for inverting said first time-dependent functional characterization into said 
30 controller; 



means for deriving a second said time-dependent functional characterization from said 
first time-dependent functional characterization, said second functional characterization 
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having a second dead-time attribute in first predefined diminishing offset from said first 
dead-time attribute, a second steady-state attribute in second predefined diminishing offset 
from said first steady-state attribute, a second curvilinear portion in first predefined 
functional offset from said first curvilinear portion, a second homaloidal portion in 
5 extrapolation of said first homaloidal portion, and a second ramped portion in 
extrapolation of said first ramped portion; 

means for deriving a third said time-dependent functional characterization from said first 
time-dependent functional characterization, said third functional characterization having a 
10 third dead-time attribute in first predefined superadditive offset from said first dead-time 
attribute, a third steady-state attribute in second predefined superadditive offset from said 
first steady-state attribute, a third curvilinear portion in second predefined functional offset 
from said first curvilinear portion, and a third homaloidal portion in truncation of said first 
homaloidal portion, and a third ramped portion in truncation of said first ramped portion; 

15 

means for determining a desired change in the value of a manipulated variable in real-time 
from said second time-dependent functional characterization, said third time-dependent 
functional characterization, the magnitude of at least one controlled variable, and said 
controller; and 

20 

means for implementing said desired change to modify said manipulated variable. 

The invention additionally provides a computer-implemented system for controlling the 
operation of a processing apparatus having at least one independently controlled 
25 manipulated variable and at least one controlled variable responsive to at least one said 
manipulated variable and further responsive to process disturbances originating 
independently of said manipulated variable, characterized by: 

means for measuring the magnitude of each controlled variable; 

30 

a model for defining an expected variation in magnitude for each controlled variable as a 
respective function of a manipulated variable disturbance instance in each manipulated 
variable; 
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means for determining, from said model, an error value in interaction of one said 
controlled variable with all said manipulated variables in achieving a manipulated variable 
modification; 

5 

means for defining a desired change in the present value of at least one manipulated 
variable from said error value; 

means for implementing said desired change to modify said manipulated variable; and 

10 

means for determining divergent response behavior in said controlling. 

Hie invention further provides a method for computer-implemented controlling of a 
processing apparatus having at least one independently controlled manipulated variable 
15 and at least one controlled variable responsive to said manipulated variable, characterized 
by the steps of: 

defining an expected variation in magnitude for each controlled variable as a respective 
function of each manipulated variable from a set of at least two models, said model set 
20 having a dynamic response inertial characteristic; and 

implementing a change defined from said model set to modify said manipulated variable 
in said processing apparatus. 

25 The invention also provides a method for computer-implemented controlling of a 

processing apparatus having at least one independently controlled manipulated variable 
and at least one controlled variable responsive to at least one said manipulated variable 
and further responsive to process disturbances originating independently of said 
manipulated variable, characterized by the steps of: 

30 

measuring the magnitude of each controlled variable; 
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defining an expected variation in magnitude for each controlled variable as a respective 
function of a manipulated variable disturbance instance in each manipulated variable from 
a set of at least two models; 

5 determining, from said set of models, an estimated modeling error value in interaction of 
one said controlled variable with all said manipulated variables in achieving a manipulated 
variable modification; 

determining an estimated process disturbance value from said controlled variable . 
10 magnitude and said estimated modeling error value; 

defining a first portion of desired change in the present value of at least one manipulated 
variable from said estimated modeling error value; 

15 defining a second portion of desired change in the present value of at least one 
manipulated variable from said estimated process disturbance value; and 

implementing said first and second desired change portions to modify said manipulated 
variable. 

20 

The invention further provides a method for computer-implemented controlling of an 
apparatus having at least one independently controlled manipulated variable and at least 
one controlled variable responsive to at least one said independently controlled 
manipulated variable, characterized by the steps of: 

25 

defining a set of consecutive discrete intervals of time in a time-dependent function; 

introducing, in each independently controlled manipulated variable, a manipulated 
variable disturbance instance of predefined magnitude, said disturbance instance 
30 prompting a response in each controlled variable; 

measuring the magnitude of each controlled variable; 
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providing a controller; 

defining, respective to each response, at least one time-dependent functional 
characterization of said controlled variable magnitude over said set of consecutive discrete 
5 intervals of time on a time-axis, said functional characterization having a zero-time time- 
axis attribute, a maximum-time time-axis attribute, a dead-time time-axis attribute, a 
response gain attribute, a ramp-rate attribute, a steady-state time-axis attribute, a 
curvilinear portion disposed between the time-axis position of said dead-time time-axis 
attribute and the time-axis position of said steady-state time-axis attribute, a ramped 

10 portion disposed between the time-axis position of said dead-time time-axis attribute and 
the time-axis position of said maximum-time time-axis attribute, and a homaloidal portion 
disposed between the time-axis position of said steady-state attribute and the time-axis 
position of said maximum-time time-axis attribute, each discrete time interval for one said 
response having the same time duration, said homaloidal portion having a zero value for 

1 5 an integrating controlled variable response, said ramped portion having a zero value for a 
non-integrating controlled variable response, and each functional characterization for one 
said response having identically valued zero-time time-axis attributes, response gain 
attributes, ramp-rate attributes, and maximum-time time-axis attributes; 

20 acquiring a first said time-dependent functional characterization respective to the response 
from measuring an effected change in said magnitude of a controlled variable after 
introducing said disturbance instance, said first functional characterization having a first 
said dead-time time-axis attribute, a first said steady-state time-axis attribute, a first said 
curvilinear portion, a first said homaloidal portion, and a first said ramped portion having 

25 its functional derivative equivalent to said ramp-rate attribute at said maximum-time time- 
axis attribute; 

inverting said first time-dependent functional characterization into said controller; 

30 deriving a second said time-dependent functional characterization from said first time- 
dependent functional characterization, said second functional characterization having a 
second dead-time attribute in first predefined diminishing offset from said first dead-time 
attribute, a second steady-state attribute in second predefined diminishing offset from said 
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first steady-state attribute, a second curvilinear portion in first predefined functional offset 
from said first curvilinear portion, a second homaloidal portion in extrapolation of said 
first homaloidal portion, and a second ramped portion in extrapolation of said first ramped 
portion; 

5 

deriving a third said time-dependent functional characterization from said first time- 
dependent functional characterization, said third functional characterization having a third 
dead-time attribute in first predefined superadditive offset from said first dead-time 
attribute, a third steady-state attribute in second predefined superadditive offset from said 
10 first steady-state attribute, a third curvilinear portion in second predefined functional offset 
from said first curvilinear portion, and a third homaloidal portion in truncation of said first 
homaloidal portion, and a third ramped portion in truncation of said first ramped portion; 

determining a desired change in the value of a manipulated variable in real-time from said 
15 second time-dependent functional characterization, said third time-dependent functional 
characterization, the magnitude of at least one controlled variable, and said controller; and 

implementing said desired change to modify said manipulated variable. 

20 The invention additionally provides a method for computer-implemented controlling of a 
processing apparatus having at least one independently controlled manipulated variable 
and at least one controlled variable responsive to at least one said manipulated variable 
and further responsive to process disturbances originating independently of said 
manipulated variable, characterized by the steps of: 

25 

measuring the magnitude of each controlled variable; 

a model for defining an expected variation in magnitude for each controlled variable as a 
respective function of a manipulated variable disturbance instance in each manipulated 
30 variable; 

determining, from said model, an error value in interaction of one said controlled variable 
with all said manipulated variables in achieving a manipulated variable modification; 
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defining a desired change in the present value of at least one manipulated variable from 
said error value; 

5 implementing said desired change to modify said manipulated variable; and 
determining divergent response behavior in said controlling. 

Further details, objects, features, and advantages of the invention are appreciated from 
10 consideration of the figures and detailed description. 

Figure 1 presents a system overview. 

Figure 2 exhibits multi-variable controller logical detail. 

15 

Figure 3 shows a basic control system block diagram. 
Figure 4 demonstrates a multi-variable control system block diagram. 
20 Figure 5 presents a time-dependent non-integrating-model characterization. 
Figure 6 presents a multiple-model non-integrating-model characterization. 
Figure 7 shows non-integrating-model set inertial characteristics. 

25 

Figure 8 exhibits a time-dependent integrating-model functional characterization. 
Figure 9 demonstrates a multiple-model integrating-model characterization. 
30 Figure 10 presents integrating-model set inertial characteristics. 

Figure 1 1 shows multi-variable control general actions in deploying a multi-model system. 
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Figure 12 presents model construction steps. 

Figure 13 exhibits dead-time time-axis attribute determination in a model. 
5 Figure 14 demonstrates Steady-state time : axis attribute determination in a model. 
Figure 15 shows ramp-portion determination detail in a model. 
Figure 16 presents curvalinear-portion determination detail in a model. 

10 

Figure 17 presents controller operation detail. 
Figure 18 exhibits adaptation methodology detail. 
15 Figure 19 shows future controlled variable requirement definition detail. 
Figure 20 demonstrates controlled variable prediction detail. 

Figure 21 presents output from a simulator for a regular DMC operating in a situation of 
20 model mismatch in modeling parameters. 

Figure 22 presents output from a simulator for the robust controller of the preferred 
embodiments operating in the situation of model mismatch of Figure 21. 

25 Figure 23 presents output from a simulator for a regular DMC operating in a situation of 
model mismatch in gains. 

Figure 24 presents output from a simulator for the robust controller of the preferred 
embodiments operating in the situation of model mismatch of Figure 23. 

30 

Figure 25 presents output from a simulator for a regular DMC operating in a situation of 
controller model inversion. 
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Figure 26 presents output from a simulator for the robust controller of the preferred 
embodiments operating in the situation of model inversion of Figure 25. 

In describing the preferred embodiments, reference is made to "logic" and to "logical 
5 portions" which are in data communication with each other. In this regard, computer- 
implemented logic generally references virtual functional elements of computer- 
implemented machine code which primarily perform tasks such as reading data, writing 
data, calculating data, performing decision operations related to data, and storing data. As 
a design consideration, the discussion of these elements as logical subportions enables a 
10 useful separation of the overall logical system into focal abstracted subcomponents which 
can each be efficiently considered, designed, studied, and enhanced within a separately 
focused and distinctively particularized context. As should be apparent, some of these 
subcomponents represent distinctive areas of specialty in their own right, even as they are 
incorporated into the comprehensive and holistic system of the described embodiments. 

15 

The specification also references the term '"Real-Time" (real-time, real time, Real-time); to 
facilitate clarity, the following paragraph presents a discussion of the Real-Time concept. 

Real-time computer processing is generically defined as a method of computer processing 
20 in which an event causes a given reaction within an actual time limit and wherein 
computer actions are specifically controlled within the context of and by external 
conditions and actual times. As an associated clarification in the realm of process control, 
real-time computer-controlled processing relates to the performance of associated process 
control logical, decision, and quantitative operations intrinsic to a process control decision 
25 program functioning to monitor and modify a controlled apparatus implementing a real- 
time process wherein the process control decision program is periodically executed with 
fairly high frequency usually having a period of between 10 ms and 2 seconds, although 
other time periods are also utilized. In the case of "advanced" control routines (such as the 
controller of the described embodiments) where a single solution instance requires more 
30 extended computational time, a larger period is essentially necessary (frequency in 

determination of changes in control element settings should be executed at a frequency 
equal-to-or-less-than the frequency of relevant variable measurement); however, an 
extended period for resolution of a particular value used in control is still determined in 
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real-time if its period of determination is repetitive on a reasonably predicable basis and is 
sufficient for utility in adaptive control of the operating mechanical assembly. 

The specification also references the term "Steady-State" (steady-state, steady state, 
5 Steady-state); to facilitate clarity, the following paragraph presents a discussion of the 
Steady-State concept. 

The solution method of the present invention also includes a general reference to the 
concept of "steady state operation". "Steady state operation" is essentially a situation 

1 0 where ( 1 ) a process is dynamically regular and uniform in its operation over a time 

interval, (2) momentum, mass, and energy entities flowing into the process are essentially 
equal to the momentum, mass, and energy entities flowing out of the process, and (3) 
accumulations of momentum, mass, and energy within the process are essentially not 
occurring unless they are explicitly expected and factored into the relevant dynamic 

15 model. Solutions of the mathematical balances with respect to the status of steady state 
operation need to also accommodate expected chemical reactions. A system in "steady 
state" is, therefore, characterized by dynamic balance, stability, steadiness, and 
equilibrium. 

20 The preferred embodiments provide a computer-implemented system for controlling the 
operation of a processing apparatus having at least one independently controlled 
Manipulated Variable (MV) and at least one Controlled Variable (CV) responsive to the 
Manipulated Variable (further responsive to process disturbances originating 
independently of the Manipulated Variable). 

25 

In defining models, circuitry and logic are provided for defining a set of consecutive 
discrete intervals of time in a time-dependent function. A "classic" DMC (Dynamic 
Matrix Control) model (a Primary Model) is a beginning point for defining an expected 
variation in magnitude for each Controlled Variable as a respective function of a 
30 Manipulated Variable disturbance instance in each Manipulated Variable. Computer 
circuitry and logic are then used in determining, from the model, an error value in 
interaction of one Controlled Variable with all the defined Manipulated Variables in 
achieving a Manipulated Variable modification. Computer circuitry and logic are also 
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provided for defining a desired change in the present value of at least one Manipulated 
Variable from the error value. Computer circuitry and logic then use the Primary Model 
(also denoted as a Reference Model) in defining a set of at least two models (a Fast Model 
and a Slow Model) for establishing an expected variation in magnitude for each Controlled 
5 Variable as a respective function of a Manipulated Variable disturbance instance in each 
Manipulated Variable; the model set (the Primary Model, the Slow Model, and the Fast 
Model collectively) also has a dynamic response inertial characteristic for each of the 
"Slow" and 'Tast" "sides" of the Primary (Reference) Model. The Reference Model, Fast 
Model, and Slow Model are also respectively referenced as the REFERENCE Model, the 
10 FAST Model, and the SLOW Model herein. 

In greater detail, the multi-variable controller of the preferred embodiments has, respective 
to each response case, at least one time-dependent functional characterization of the 
Controlled Variable magnitude over the set of consecutive discrete intervals of time on a 

15 time-axis, the functional characterization having a zero-time time-axis attribute, a 
maximum-time time-axis attribute, a dead-time time-axis attribute, a response gain 
attribute, a ramp-rate attribute, a steady-state time-axis attribute, a curvilinear portion 
disposed between the time-axis position of the dead-time time-axis attribute and the time- 
axis position of the steady-state time-axis attribute, a ramped portion disposed between the 

20 time-axis position of the dead-time time-axis attribute and the time-axis position of the 
maximum-time time-axis attribute, and a homaloidal portion disposed between the time- 
axis position of the steady-state attribute and the time-axis position of the maximum-time 
time-axis attribute, each discrete time interval for one response case having the same time 
duration, the homaloidal portion having a zero value for an integrating Controlled Variable 

25 response, the ramped portion having a zero value for a non-integrating Controlled Variable 
response, and each functional characterization for one response case having identically 
valued zero-time time-axis attributes, response gain attributes, ramp-rate attributes, and 
maximum-time time-axis attributes. 

30 The first time-dependent functional characterization respective to the response (the 
Primary Model) is derived from measuring an effected change in the magnitude of a 
Controlled Variable after introducing the disturbance instance. The first functional 
characterization has a first dead-time time-axis attribute, a first steady-state time-axis 
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attribute, a first curvilinear portion, a first homaloidal portion, and a first ramped portion 
having its functional derivative equivalent to the ramp-rate attribute at the maximum-time 
time-axis attribute. 

5 The second time-dependent functional characterization (the Fast Model) is derived from 
the first time-dependent functional characterization, the second functional characterization 
having a second dead-time attribute in first predefined diminishing offset from the first 
dead-time attribute, a second steady-state attribute in second predefined diminishing offset 
from the first steady-state attribute, a second curvilinear portion in first predefined 
10 functional offset from the first curvilinear portion, a second homaloidal portion in 

extrapolation of the first homaloidal portion, and a second ramped portion in extrapolation 
of the first ramped portion. 

The third time-dependent functional characterization (the Slow Model) is also derived 
1 5 from the first time-dependent functional characterization, the third functional 

characterization having a third dead-time attribute in first predefined super-additive offset 
from the first dead-time attribute, a third steady-state attribute in second predefined super- 
additive offset from the first steady-state attribute, a third curvilinear portion in second 
predefined functional offset from the first curvilinear portion, a third homaloidal portion in 
20 truncation of the first homaloidal portion, and a third ramped portion in truncation of the 
first ramped portion. 

The Controller for the multi- variable controller is defined by essentially inverting the first 
time-dependent functional characterization (the Primary Model). 

25 

In implementing the models, an inertial characteristic value is input into the database of 
the model-variable controller so that the dynamic response inertial characteristics are 
achieved in the model set. The inertial characteristics (establishing robustness in 
response) are (a) different for the "Slow side" of the Primary Model and for the "Fast 
30 side" of the Primary Model, or (b) the two inertial characteristics are identical. In this 

regard, the first, second, and third time-dependent functional characterizations define three 
models in a model set, the model set incorporating a first dynamic response inertial 
characteristic between the first and second time-dependent functional characterizations 
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and a second dynamic response inertial characteristic between the first and third time- 
dependent functional characterizations. 

After the inertial characteristics have been input, the Primary Models are defined by 
5 introducing, in each independently controlled Manipulated Variable, a Manipulated 
Variable disturbance instance of predefined magnitude, the disturbance instance 
prompting a response in each Controlled Variable. The magnitude of each Controlled 
Variable is then measured, retained, and used in definition of the Primary Model according 
to traditional DMC practice (reference previously discussed US Patent 4,349,869). 

10 

In implementing use of the multivariable controller, a number of different approaches are 
possible. 

In one approach, a change is defined from the model set to modify each Manipulated 
15 Variable in the processing apparatus; in achieving this definition, a desired change is 
determined for the value of a Manipulated Variable in real-time from the second time- 
dependent functional characterization, the third time-dependent functional 
characterization, and the magnitude of at least one Controlled Variable. 

20 Alternatively, a desired change is determined for the value of a Manipulated Variable in 
real-time from the first time-dependent functional characterization, the second time- 
dependent functional characterization, the third time-dependent functional 
characterization, and the magnitude of at least one Controlled Variable. 

25 In a further embodiment, an estimated modeling error value in interaction of one 
Controlled Variable with all the Manipulated Variables is determined in achieving a 
Manipulated Variable modification. An estimated process disturbance value is then 
determined from the Controlled Variable magnitude and the estimated modeling error 
value; and a first portion of desired change in the present value of at least one Manipulated 

30 Variable is determined from the estimated modeling error value. Then a second portion of 
desired change in the present value of at least one Manipulated Variable is determined 
from the estimated process disturbance value. The first and second desired change 
portions are finally effectively combined to modify each Manipulated Variable. 
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A basis for adaptation of the models in real-time operation is also enabled in the described 
embodiments. In this regard, the modeling error value and process disturbance values are 
used in decisions to adapt the controller in the context of the inertial characteristics of the 
5 model set. Furthermore, determination and counteraction respective to divergent response 
behavior in the controlling operation are implemented to provide strategic robustness via 
(a) the adaptive process and (b) ongoing use of the models which have been adaptively 
modified in the controller. 

10 A formalized characterization of the multi- variable controller is set forth in the following 
section. The discussion will be divided into six subsections: a first subsection describing 
the models of the robust multi-variable controller; a second subsection discussing the 
controller definition in the robust multi- variable controller; a third subsection discussing 
determination of the process in the robust multi-variable controller; a fourth subsection 

15 discussing the adaptive process in the robust multi- variable controller; a fifth subsection 
discussing the decision to adapt in the robust multi- variable controller; and a sixth 
subsection discussing determination of divergence and adaptation to divergence in the 
robust multi- variable controller. 

20 The following first subsection of the formalized characterization discusses the models of 
the robust multi- variable controller. 

The robust matrix process control method of the preferred embodiments requires sets of 
numerical values, which explicitly define the relationship between (a) input variables, 

25 including feedforward variables (denoted as manipulated variables, MV), and (b) output 
variables (denoted as controlled variables, CV). Each data set is referred to as a "model". 
Each data set is expressed as a vector in which the entries are the values of the data set in 
ordered relation to time. Note that the input variables to the model (MV) are the output 
signals from the control system (MV) and that output variables from the model (CV) are 

30 input signals to the control system (CV). 

These data sets (expressed as vectors) are processed by computer-implemented procedures 
to (a) produce the mathematical datalogical objects (scalars, vectors, and matrices) and (b) 
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to achieve engine-solved computations (implementation and manipulation in computer 
logic of sets of scalars, vectors, matrices, and mathematical operations which are arranged 
in such an order as to perform a non-trivial mathematical calculation) to achieve robust 
control of the controlled apparatus. 

5 

As noted previously, the term "model" references a set of time stamped values defined as 
a vector 

Equation 1 

10 

M(t)= {a (t) }, where t=0, l*At, 2*At, 3*A,...np*At 

where "np" is the model size - therefore, there are "np" a( t ) values. 

15 Each value a< t ) is obtained through regular model identification techniques. This usually 
involves plant testing and data collection techniques along with subsequent analysis and 
processing. M(t) is alternatively a step model, an impulse model, or a parametric model. 
The description below refers to step models, but an equivalent exists for either an impulse 
model (an impulse model is readily transformed into a step model by a summation 

20 process) or a parametric model (transformed to either a step model or impulse model). 

There is always a model M(t) for each pair of M V and CV. The trivial model M(t)= {0} 
from t=0 to np*At, indicates that no physical relation exists between the two elements of 
the pair, but this model exists nonetheless. Usually the model is represented after 
25 smoothing of the sampled data (that is, each a( t) ) to obtain a continuous curve; the 

smoothed sampled data is shown in Figure 5 and discussed in the specification respective 
to Figure 5. 

A set of models is used to define the matrix A, where each column of the matrix refers to a 
30 specific predicted future CV value. It should be noted that the "A" matrix is apparent to 
those of skill in Dynamic Matrix Control, and is normally used in descriptions of the 
technology. 
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This first set of models obtained by plant testing is denoted herein as the "Reference 
Model". Since the "Reference Model' 5 is time dependent, the time ("t") is defined as the 
"Reference Time". "Reference Time" goes from zero to t fina i. The definition of a 
Reference Time is important in deriving two other models from the Reference Model. In 
5 this regard, the preferred embodiments are based in part upon defining a set of consecutive 
discrete intervals of time in related time-dependent models where each discrete time 
interval for each related model has the same time duration and each related model also has 
identically valued zero- time time-axis attributes and maximum-time time-axis attributes. 
The Reference Model and the two additional models (further described herein) therefore 
10 all have a time axis with a time span from t=0 to the same tf ina i. A fundamental idea in the 
preferred embodiments is to define and use an "envelope" derived-from and respective-to 
the "Reference Model". The "envelope" is defined such that it delimits a zone of 
uncertainty where the real model is trusted (a) to exist and, further, (b) to evolve with time 
and process conditions. 

15 

Thus, the "envelope" has an "upper" boundary and a "lower" boundary (respective to the 
time axis); these boundaries require definition. For simplicity, assume that each model 
M(t) has a continuous mathematical equivalent. That is: M(t) is a z-transform of M*(t)(M* 
being continuous). Thus: 

20 

Equation 2 

M(t)= Z(M*(t))= aO + al*z _1 + a2*z 2 + etc.. 

25 Each model exhibits dead time. Therefore, M*(t) is written as M*(t-d). M* is zero when 
t<d. 

Each model has a final steady-state gain or a final constant ramp rate; these are denoted as 
G or R, respectively. The time of steady-state (or final constant ramp rate) initiation is 
30 defined as t^. The last value of the model occurs at time W 

The envelope boundaries are the lower limit model and the upper limit model for a 
positive gain model. The lower limit model is denoted as the "Slow Model" and, 
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respective to the Reference Model, has the physical significance of a process (a) having 
more dead time and (b) reaching steady state "later''. The upper limit model is denoted as 
the "Fast Model" and, respective to the Reference Model, has the physical significance of 
a process (a) having less dead time and (b) reaching steady state "sooner". The 
5 fundamental idea of this robust controller is that, as long as the controlled process exhibits 
a behavior that "lies" between the "Slow Model" and the l Tast Model", then it is possible 
for a controller to behave like a perfect controller in the sense that there will be no model 
mismatch for which a compensation must be calculated and implemented. 

10 In defining the "Fast Model" and the "Slow Model", substitutions of the time variable into 
M*(t) are performed; the aim of the substitution is to contract (Fast Model) or expand 
(Slow Model) the function of the Reference Model along the time axis. 

The "Fast Model" is obtained by the following substitution of the time variable into M*(t): 

15 

For time before Dead Time: 
Equation 3 

20 t = t fast /(L0-F rl ) yielding t fast = (1.0-F rl )*t . 
For time after Dead Time: 
Equation 4 

25 

tfast = (1.0.F rl )*t d + (1.0-F r2 )*(t-t d ) . 

For a ramp model, F r2 is zero. 

30 The portion after tfast>tss is filled with the model value at steady-state up to tfimi. In the 
case of a ramp, the ramp is extended at the constant rate R. 
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The "Slow Model" is obtained by the following substitution of the time variable into 
M*(t): 

For time before Dead Time: 

5 

Equation 5 

tsIow=(1.0+F r3 )*t . 

10 For time after Dead Time : 
Equation 6 

tsbw = (1.0+F r3 )*t d + (L0+F r4 )*(t-td) . 

15 

For a ramp model, is zero. 

The portion after t s i ow >tf ma i is dropped. 

20 F r i, Fr2, F r 3, and F r4 are robustness factors. The robustness factors are positive adjustable 
values usually in the numeric range of 0-.5, but they can be of greater magnitude. The fact 
that the time-referenced model contraction and expansion are treated differently before 
and after dead time allows for more flexibility in building the limit models (the "Slow 
Model" and the "Fast Model" are collectively "limit models" defining the "envelope" 

25 around the "Reference Model"). In this regard, (a) the veritable dead-time of the 
operational apparatus being controlled and (b) the veritable transient trajectory of the 
operational apparatus being controlled (as characterized in each model with a "curvilinear 
portion" as further described herein) both are subject to time-variance and/or to 
measurement error; the flexibility afforded by the independent methods (1) in "Reference 

30 Model" expansion to the "Slow Model" and (2) in "Reference Model" contraction to the 
"Fast Model" facilitates the robust operational response enabled by the preferred 
embodiments. Furthermore, the "Slow" and "Fast" Models need this flexibility when 
either of these limit models solves to be essentially proximate to the Reference Model. 
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At the end of the above substitution processes, all time in the Reference and limit models 
will have the same equivalence, with each model providing a time axis of a set of 
consecutive discrete intervals of time where each discrete time interval for each related 
5 model has the same time duration and each related model also has identically valued zero- 
time time-axis attributes and maximum-time time-axis attributes, as previously noted. 
That is: 

Equation 7 

10 

t = tfest = tslow • 

As should be appreciated, before processing all the models in the dynamic control matrix 
to produce all the affiliated envelopes, some or all of the Reference Models are optionally 
15 extended with a higher magnitude of W to ensure that the last operation keeps unchanged 
(a) the gain of any non-integrating process model or (b) the ramp rate of any integrating 
process model. 

The values of F r i, F^, F r 3, and F r 4 belong to a model and are intra-functional. It should be 
20 noted that it is possible (even desirable) to have fewer values to adjust than represented by 
the full number of robust factors. In this regard, note that the mathematics herein 
described allow for any set of F r factors to be grouped and assigned a single value through 
use of ratios and user-defined functions. These F r values are local to each controlled 
variable and/or each manipulated variable; it is therefore possible to provide a single F r for 
25 adjustment respective to the whole controller. This feature does not prevent local 
adjustment if desired; therefore, each F^ applied to each model is decomposed as: 

Equation 8 

30 Fric(i j)= F^global) + F r (whole MVj) + F r (whole CVj) + Frktspecific to a Reference Model 
instance) . 
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The requirement that each F^i j) be a positive number is still valid, but this does not 
prevent other components in the above equation from having a negative value. In this 
regard, a local F r may be negative with the effect that the net Frk(i j) is of lower magnitude 
than FXglobal) even as Frfc(ij) is positive (or even zero). 

5 

When all instances of F r (i j) have been chosen, the models have been properly extended (if 
necessary), and M*^ and M* s i ow have been calculated, then M*fest and M*siow are written 
into sampled form (or Z-transform form). The diagram and discussion respective to 
Figure 16 illustrates a computer-implemented approximative process which yields 
10 equivalent results to the above formalized models. The process involves interpolation of 
the Reference Model M (already in z-transform form). 

In vector form an example result is: 

15 Reference Model (0,0,0,0,0,0,ai,a2,a3,a4,a5, ... a^, ... a^) 
Slow Model (0,0,0,0,0,0,0,bi,b2,b 3 ,b4,b5, ...ass, ... a ss ) 
Fast Model (0,0,0,0,ci,C2,c 3 ,C4 5 C5 5 . . . a^, . . . a ss ) 
where b\ and Ci are interpolated from ai with the rule: 

20 Equation 9 

Xi = f*ai + (1.0-f)*ai + i 

and f is calculated from each F* and current time t in the models. 

25 

It should be emphasized that there are user choices in defining and calculating limit 
models because these models are directed to defining the boundary (limits) of controlled 
process response (behavior) where the boundary separates (a) an operational zone of 
predicted behavior where the model is trusted to be exact from (b) a complimentary 
30 operational zone where the model is trusted to have detected a disturbance in the process. 
Therefore, a benefit of the described embodiments is realized insofar as imprecision in the 
Reference Model does not diminish robustness to the degree that similar imprecision in the 
Reference Model diminishes robustness when there is not a defined envelope around the 
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Reference Model (that is, as in traditional DMC). Also noteworthy is the fact that the 
proposed transformation touches only (a) the dead time and (b) the process model 
trajectory; the gains of the model are unchanged. Adaptation for the gain is added later; 
but adaptation is not required to achieve robustness in the face of (a) model imprecision, 
5 (b) (more importantly) dead time uncertainties, and/or (c) dead time variations respective 
to either (1) time or (2) process conditions. 

The following second subsection of the formalized characterization discusses the 
controller definition in the robust multi-variable controller. 

10 

The overall structure of the controller proper is the same as for a regular DMC (Dynamic 
Matrix Controller). In this regard, the controller proper has three blocks: the predictor 
block, the optimizer block, and the move calculator block. However, as will be further 
discussed in the discussion of the preferred embodiments, the predictor and the move 
15 calculator blocks differ from traditionally deployed DMC. The optimizer block is a 
conventional LP solver or any other commercially available steady-state optimizer. 

The predictor block of the preferred embodiments differs from traditional DMC in using 
the Slow and Fast Models with the optional use of the Reference Model. In this regard, 
20 the important element of difference relates to the prediction of at least two limit behaviors 
of the process. 

The controller block uses the usual least-squares equation of traditional DMC to evaluate 
the final moves in the manipulated variables; but, in the preferred embodiments, the 
25 calculation of E(i) is different from the previous art in DMC - in this regard, the traditional 
calculation is: 

Equation 10 

30 AU=(A t WA) _1 .A t W^E 

where AU = MV t+1 - MV t . 
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(A'WA)' 1 is also denoted as the "ATA" matrix herein and A'WE is also denoted as the 
"ATE" matrix herein. 

As in the case of traditional DMC, "move suppression factor", "equal concern error", and 
5 "constraints handling" are calculated and processed in the preferred embodiments when 
solving this least-squares equation. These are independent of the vector E. In traditional 
DMC, robustness is substantively achieved with "move suppression commissioning" 
under the assumption that the best possible (Reference) model was obtained. This is 
implemented, however, at the expense of suppressing the reactive movement of the 
10 manipulated variables. Consequently, the controller may have diminished responsiveness 
to disturbances and set point changes respective to fully desired responsiveness. 

But controller robustness can be obtained as long as the controlled process exhibits the 
predicted process behavior. Since exact prediction of the controlled process is impossible, 
15 the determination of the process is the location where robustness is "found". So, for a 
controller to behave like a perfect controller in the sense that there will be no model 
mismatch to compensate-for, the model mismatch has to be minimized by a modification 
of the process determination. The above control equation is therefore modified in the 
preferred embodiments to be: 

20 

Equation 1 1 

AU= (A t WA)' l -A t W.F(models l Y,-.) 

25 where F is a function dependant of process measurements Y, and errors based on limits 
models are used. In this regard, F(models,Y,..) has replaced E in the 
AU= (A'WA)" 1 - A l W-E equation. 

The following third subsection of the formalized characterization discusses determination 
30 of the process for the robust multi- variable controller. 

For each CV there are at least two predictions: the one from the Fast Model CV^ and the 
one from the Slow Model CV s i 0W . Optionally, the prediction from the Reference Model is 
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used along with at least one limit model. The predictions are obtained as in traditional 
DMC practice with: 

Equation 12 

5 

CV(i, = (oft j, l) * HUU, k + / - 1) ■ 

However, in the preferred embodiments, there are at least two prediction blocks: one for 
the Slow Model and one for the Fast Model. It is also preferable to maintain the 
10 prediction from the Reference Model, but, again, this is optional. 

The measured output Y is compared to these predictions through an ad-hoc function. The . 
function is designed to generate a modified measured output Y*. Then the "error to 
correct 55 is a function of at least these three elements and, optionally, the Reference Model 
15 prediction: 

Equation 13 

Y= f(CV fest (t),CV slow (t),Y,etc.) . 

20 

In example, the function is defined as: 
Equations 14-16 

CV mx = Max(CV fast (t-k), CVs, ow (t-k) ) {Vk | -d<k<+d} 
25 CV min = Min(CV f ast(t-k), CV s!ow )(t-k) ) d= 0,1,2 or 3 (chosen) 
CV mean = (CV ftlt +CV ltow )/2.0 

or alternatively: 

30 
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Equation 17 



CV me an= CV re f erC nce at the expense of more computing, but it might be preferred insofar 
as a better estimate is achieved. 

5 

The extension { Vk | -d<k<+d} within the predicted values expands the choice of 
predicting models to those having less and more dead time. So, if d=3, for example, then 
this function actually uses at least 14 different models, with each model generating 
different responsive behavior: 

10 

Equations 18-20 



If Y>CV fast then Y*=Y-CVfast +CVfast+CV sl0W 
If Y<CVfcs t then Y*=Y- CV m +CV^ 0V/ 
15 Otherwise Y*= CV mca n • 

Another possible function (that has the benefit of being continuous) is the following as 
obtained by defining a function <S> to normalize variations in the predicted CV's: 

20 Equation 21 

<D = — I 2 _ ; 

it is possible to calculate the desired Y* according to 

25 

Equation 22 



Y =- 



cv fn „ + cv. 



slow 



fast 





f 




* 




V 



/ l+a«D") 
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10 



The value of "a" and "n" are adjusted to provide a behavior that fits statistical expectation. 
For example, a=l .5 and n=4 produces a function that has desirable effects: 

it returns Y* equal to the CV mean when Y is within CV^st and CV s i 0W ; 

it returns Y* equal to the measured CV when far away from CVfest or CV sIow ; and 

it returns Y* equal to the measured CV mea n plus a small correction when close to CVfest or 

CV s iow 

The function term 



Equation 23 



f«-— i 

^ l + a®" ) 



15 is a confidence function that is zero (or close to zero) when the measured output is within 
the Slow and Fast predictions. In the first example, the confidence interval is defined by 
CV max and CVmin. In both cases, the function is a scaling sensitivity factor. The function 
also jumps (increases) toward unity when the measured output moves away from the 
predictions; this indicates a confidence that the measured output indicates a strong 

20 deviation in the plant (apparatus in operation) needing correction. Also, CV mean (i.e, (CVfest 
+ CVsiow )/2 ) is optionally replaced by CV re ference (if available). 

The error to compensate-for is now estimated from Y*- CV^^ for each CV. This 
formulation now contains the process information from the Slow and Fast Models, and it 
25 generates a controller using at least 2 models simultaneously (and, optionally, more 
models as previously indicated): 

Equation 24 



30 AU= (A'WA)" 1 - A'WCY'-CVn^) 
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The following fourth subsection of the formalized characterization discusses the adaptive 
process in the robust multi-variable controller. 

5 Further robustness can be added to the controller. Since the robust controller is 

presumably stable "in the face of (when dealing with or when subject to coping with) 
process variation caused by dead time change and overall model change, then the 
(remaining) gain element is now adapted on-line (during real-time operational use of the 
preferred embodiments) with better efficiency respective to traditional DMC practice since 

10 the overall control process of the described embodiments presumably provides relatively 
improved stability. This adaptive process is perhaps more precisely denoted as "semi- 
adaptive" insofar as (a) only the gain is adapted and also (b) because some discrepancy is 
inherently accepted in model shape and model dead time by the robust controller structure 
described by Equation 1 1 . The information derived from using two limit predictions is 

15 further used for estimating (on-line) improved gains for all the models (Reference, Slow, 
and Fast Models). By guiding the adaptation process to operational zones where 
maximum stable information is available, the adaptive algorithm is fed with reliable and 
non-disruptive data. This enables the gain estimation algorithm to proceed smoothly. 

20 The process of gain adaptation involves comparing past prediction of the Reference Model 
with historical data from the process. Then multiplicative correction factors are estimated 
by the least-squares method. 

The model prediction based on the Reference Model at time t k for each CV(i) can be 
25 written as: 

Equation 25 



and this is transformed into a larger sampled time interval and into a difference equation: 




30 
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Equation 26 

/-I M 

5 where: N is the number of MV's and FF's 

"np" being the model size; then there are "r" divisions of equal size "s" in the model such 
that there are np = r*s values in the model. 

10 ACV s (i,k) and AU s (i,k) are taken as the time difference of their values separated by a time 
interval equal to s*At. Usually, "s" is chosen such that effect of dead time is minimized; 
this generates stability in the adaptation process. 

The ACV s (i,k) values are weighted according to their age and reliability. Age is 
15 discounted by weighting ACV s (i,k) with X(t) so that recent values have more weight in the 
adaptation calculation than relatively older data: 

Equation 27 

20 X w =^*- |) 

where Xo is adjusted to lie usually between .9 and 1 .0. 

The reliability is related to the model predictions and, therefore, involves the predicted 
25 values obtained by the Slow and Fast Models as well as the type of MV movements. It 
also uses on-line process value Y. Since a Fast Model and a Slow Model each define a 
zone of uncertainty respective to the Reference Model, this overall information is used to 
determine readiness of the data for adaptation use: 

30 Equation 28 



V (W =f(CV sl ow, CVfest, CV reference , MV, Y) and 0.0 < V (i>k) <1.0, 
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that is,: 

Equations 29-33 

5 

IfCVslow^CVfest 

if Y>CVsiowandY>CVfcst then V(i,k)=0.0, 
if Y<CV sl0W andY<CV fest then V(i,k)=0.0, 

10 otherwise 

V(i,k) = MinCIY-CV^UY-CVsUyiCVfasrCV^w)!, 

otherwise 

15 

V(i,k)=1.0 . 

Finally, if AU s G',k) is not monotonic within its "s" time interval according to its historical 
record, then V(i,k)= 0.0. 

20 

The monotonicity of U is evaluated by: 

Equation 34 

25 |AUs|<Ski|AU t+ i| 

with Zk» = L0 with each kj > 0 (chosen to preferentially weight the effect of first AU t +i 
values in the historical data series). 

30 Optionally, a part of the proposed scaling sensitivity factor determines the V function: 
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5 Whatever the choice, this function is chosen by the user to eliminate (or minimize) the 
effect of model mismatch caused by model parameters other than gain by using the 
predictions obtained from at least two limit models. The selected function varies according 
to the characteristics of the process such as (a) a deterministic process with few 
disturbances or (b) a process characterized by strong stochastic disturbances. The function 
10 V has the effect of (a) screening data containing the most valuable information for 
adapting the process gains and (b) rejecting some of the transients and disturbances 
introducing errors in the adaptation process. 

Continuing, each ACV s (i X) is modified by the weighting equation: 

15 

Equation 36 

ACV s (i,k)= ACV s (i,k)*( X {k) * V (U0 ) . 

20 This forms the matrix " C " used to solve the following minimization of sum-of-squared 
error: 

Equation 37 



25 s(o=IX, *v(i,k))*E{w 2 =2>(*)**(a) 2 > 



the summation being done on available historical data 



30 
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Equations 38-39 

E(U k) = ft k) - ACV S ft *) = 
Ar,ft*)-££AG(UW^^ 

y=i M 

5 To minimize S(i), a set of changes in gain is introduced to provide multipliers of the actual 
"a" values of each model. These are the AG(i j) multiplying each term of the summation 
(Equation 39). Therefore, there are as many AG as there are models. For future reference, 
AU and ACV quantities are collectively denoted as "delta-s quantities". 

10 Continuing, the changes in gain are computed from the least-squares equation: 

Equation 40 

AG(i)= (C'HCy'-C'HAYfl) (repeated for each CVj). 

15 

The resulting individual AG(j) multipliers are validated and relaxed (filtered) so that 
smooth changes are implemented. In some cases, the validation process results in the 
rejection of the entire vector estimate. Note that, in this case, the validation process adds 
robustness to the adaptation. 

20 

Continuing, each Model (Slow,* Fast and Reference) is updated by the matrix AG(i): 
Equation 41 
25 M= M AG(i) . 

This is repeated for each selected CV(i). Each CV(i) does not necessarily need adaptation 
in this regard; this is a choice (a) defined at controller commissioning or (b) implemented 
in real-time operation by being either "turned on" or "turned off" respective to each 
30 individual CV. 
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/The following fifth subsection of the formalized characterization discusses the decision to 
y adapt in the robust multi-variable controller. 

The implementation of the results of adaptation: 

5 

Equation 42 

AG = (C t HCy 1 .C t HAY 

10 affects all models and may produce either an ill-conditioned matrix or a singular matrix. 
As should be apparent, a singular matrix or an ill-conditioned matrix should not be 
implemented for control. Therefore, the gains are validated by analyzing the 
consequences of any gain change. This adds further robustness into the adaptation 
process. 

15 

The initial model is used to generate a set of characteristic values such as eigenvalues or 
singular values. The calculation is usually performed with a commercial package. 
Therefore, only the resulting vector and matrices from the decomposition calculation 
process are important here. 

20 

These characteristic values will be used to validate future changes and they serve as 
reference values. 

The following assumes that the gains have been normalized to compensate for scale and 
25 physical unit effects. The gain matrix is made a-dimensional. 

Since the matrix of the process gain is the target of the adaptation, the characteristic values 
are obtained from G r ^ erCRCt : 

30 Equation 43 

For eigenvalues: G re fererenee= T l -X-T ; 
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Equation 44 

For singular values: G re fererence = U-S-V . 

5 The matrix X or £ contains the characteristic values to monitor. For an ill-conditioned 
matrix, it is necessary (a) to detect a change of sign in one characteristic value of the 
process and then (b) to responsively implement the equivalent change of sign in the model 
gains to produce the same sign in all the characteristic values and add to the controller 
stability. 

10 ' 

The vectors T and/or U,V are maintained as constant, and only the characteristic values 
are updated. Then, as the adaptation calculates changes in model gains G(i,j), new 
corresponding X or E are determined and validated to produced A,* or E*. The validated 
characteristic values are then used in back-calculating the final gains: 

15 

Equation 45 

For eigenvalues: G* rc f jS io W> fast = T X X*-T ; 
20 Equation 46 

For singular values: G* re f^iow,fast = U-Z**V . 

It should be noted that some safeguards are preferably introduced to prevent abrupt 
25 changes. In this regard, large characteristic values (greater than or equal to threshold e) 
and small characteristic values (less than threshold e) are handled in different ways. The 
small values are associated with instabilities and potential divergence because the 
likelihood of having sign difference between the controlled process and the model is 
greater. The small characteristic values are accordingly not allowed to decrease in 
30 magnitude. The threshold is defined (denoted e) to prevent any magnitude change from 
crossing a limit where the change can induce divergence in the controller (that is, inverted 
gain). The divergence problem is not solved by the adaptation mechanism. A different 
mechanism is described in subsection 6. 
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If any characteristic value exhibits a sign change and (at the same time) is a large value, 
then the adaptation for the gain related to that CV will be rejected as unreliable. 
Otherwise, the characteristic values are validated and relaxed according to: 

5 

Equations 47 and 48 

X\ t = max(s,(1.0-h)*A Te ferencc > min(e, (l.O+hJ^XreferenccA*^)) 

XV XV f*(X\ t - X\t-i) with f chosen 0.0<f<1.0 

10 

where "h" defines clamping limits respect to the initial characteristic values and "f limits 
speed of changes. Optionally, "f" is adjusted based on historical data: 

Equation 49 

15 

f= g( statistics of past X) . 

The unfiltered characteristic values are maintained in history (a) for future use when "f ' 
depends on the statistics of past calculated characteristic values or (b) to assess the 
20 stability of the adaptation process. 

The following sixth subsection of the formalized characterization discusses determination 
of divergence and adaptation to divergence in the.robust multi-variable controller. 

25 As mentioned before, the initial Reference Model is used to calculate characteristic values 
such as eigenvalues and eigenvectors, or the singular value decomposition E. 

The calculation is performed with a commercial package such as "Recipes in C" by 
Cambridge University Press. 

30 

The process of adaptation to divergence does not depend on the semi-adaptative 
mechanism of the prior two subsections and runs independently. 
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The set of eigenvectors or the vectors of V (in the singular value decomposition) 
constitutes a mapping of the current MV's into a space where effects are decoupled. 

The vectors T and/or U,V (kept constant) are used to transform the ACV(i) and the 
5 AMV (i). The evaluation of these last two quantities is performed by using the "delta-s 
quantities": 

Equation 50 

10 ACV s (i,k)andAMV s (i,k) 

(once again) to minimize effect of dead time. This produces stability in the divergence 
detection mechanism. 

15 Then each AMV s (i,k) and ACV s (i,k) is validated by the validity factor V(i,k). This uses 
the information coming from the Slow and Fast Models indicating the transient nature of 
the changes in the manipulated variables and controlled variables: 

Equations 51 and 52 

20 

AMV s (i,k)=AMV s (i,k)*V (W ; 
ACV s (i,k)= ACV s (i,k)*V (iJ0 . 

Finally a transformation is done by: 

25 

Equation 53 

For eigenvalues: A=T-ACV S and B=T-AMV S ; 
Equation 54 

30 

For singular values: A=U -1 ' ACV S and B=V AMV S . 
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The elements of A are the transformed CV and the elements of B are the transformed MV. 
This introduces a decoupling of effects. 

Continuing, an estimate of the sign of the characteristic values is obtained by a ratio of the 
5 non-zero elements and/or sufficiently large elements of B to the elements of A. A "too- 
small" value of any element of A indicates, obviously, that the corresponding MVs (a) did 
not move or (b) did not move in concert to act on the corresponding transformed CV. In 
that case, no conclusion is reached regarding the corresponding characteristic value and, 
accordingly, nothing is changed. 

10 

The known small characteristic values are tracked in the historian, and a persistent sign 
change triggers a positive divergence detection signal. The logic here is primarily based 
on a user chosen algorithm that uses statistics such as those derived from a CUSUM 
algorithm. 

15 

When divergence is detected, changes in the model gains G(ij) (corresponding to the 
value of X or £, but with the sign inverted) are introduced to produced X d or E d . This 
means that (for eigenvalues) the following transformation is accomplished for any small 
eigenvalues (other eigenvalues may have been adapted as well): 

20 

Equation 55 
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25 The new vector, denoted as X with inverted characteristic values, is then used to back- 
calculate the final gains: 
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Equation 56 

For eigenvalues: G d K r fi \ aw ^ st = T l 'X d -T ; 
5 Equation 57 

For singular values: G d re fri 0 w > &st~ U*E d V . 

The new gains are then inserted into the Reference Model and any other model previously 
10 derived from it (for example, the Fast and Slow Models). It is important to maintain the 
same gains for all corresponding parts of all models. 

Further details in the preferred embodiments are appreciated from consideration of the. 
figures and their description. Turning now to the figures and Figure 1, System Overview 
15 100 presents an overview of the physical components in an applied described embodiment. 

Control Computer 104 incorporates Control Computer CPU 122 for execution of Control 
Computer Logic 120 in real-time operational monitoring and control of Operated 
Apparatus 102. Operated Apparatus 102 is alternatively, in example and without 
20 limitation, a steam turbine, a gas turbine, a chemical process, an internal combustion 
engine, or a furnace. 

Communication Interface 106 incorporates Communication Interface CPU 1 18 for 
execution of Communication Interface Logic 1 16 in facilitating bilateral data 
25 communication between Control Computer 1 04 and Multi-Variable Controller 1 08. 

Multi- Variable Controller 108 incorporates Multi-Variable Controller CPU 1 14 for 
execution of Multi- Variable Controller Logic 1 12 in the methodology of the preferred 
embodiments; Multi- Variable Controller Logical Detail 200 of Figure 2 shows further 
30 detail in Multi- Variable Controller Logic 1 12. Monitor and Keyboard 1 10 provides a 
monitor for human viewing of information and data. 
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Analog Input Signal 124 provides an analog input signal from Operated Apparatus 102 as 
an electrical voltage input to Control Computer 104 and is responsive to an attribute (for 
example, a temperature or a pressure) within of Operated Apparatus 102. After 
conversion into a datalogical entity within Control Computer 104, Analog Input Signal 
5 124 is denoted as a Controlled Variable (CV) insofar as the control of measured attributes 
of Real-time Interface Logic 210 is the goal of a controller deployed in either Control 
Computer 104 or Multi- Variable Controller 108. 

Digital Input Signal 126 provides an digital input signal from Operated Apparatus 102 as 
10 an electrical voltage input to Control Computer 104 and is responsive to an attribute (for 
example, a valve being either open or not open) within Operated Apparatus 102. After 
conversion into a datalogical entity within Control Computer 104, Digital Input Signal 126 
is denoted as a Controlled Variable (CV) insofar as the control of measured attributes of 
Real-time Interface Logic 210 is the goal of a controller deployed in either Control 
15 Computer 104 or Multi-Variable Controller 108. 

Digital Output Signal 128 provides an digital output signal from Control Computer 104 as 
an electrical voltage input to modify and manipulate the position of a control element in 
Operated Apparatus 102. In so doing, Digital Output Signal 128 modifies attributes (for 
20 example, the kinetic energy available for transfer to a fluid in enabling fluid transport 
through a pump) of Operated Apparatus 102. Accordingly, a datalogical entity within 
Control Computer Logic 120, Communication Interface Logic 116, and Multi-Variable 
Controller Logic 112 whose value is represented in the voltage of Digital Output Signal 
128 is denoted as a Manipulated Variable (MV). 

25 

Analog Output Signal 130 provides an analog output signal from Control Computer 104 as 
an electrical voltage input to modify and manipulate the position of a control element in 
Operated Apparatus 102. In so doing, Analog Output Signal 130 modifies attributes (that 
is, the cross sectional area available for fluid flow within a pipe) of Operated Apparatus 
30 102. Accordingly, a datalogical entity within Control Computer Logic 120, 

Communication Interface Logic 1 16, and Multi- Variable Controller Logic 1 12 whose 
value is represented in the voltage of Analog Output Signal 130 is denoted as a 
Manipulated Variable (MV). 
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The natures of Operated Apparatus 102, Control Computer 104, and Communication 
Interface 106 are presumed to be apparent to those of skill and are presented here to enable 
an overall understanding of the preferred described embodiments and their use. Details in 
5 Multi- Variable Controller 108 are focal in most subsequent discussion in this 

specification. An example of Control Computer 104 is given in U.S. Patent Serial No. 
09/482,386, filed January 12, 2000, entitled "PROCESS CONTROL SYSTEM WITH 
INTEGRATED SAFETY CONTROL SYSTEM". 

10 Turning now to Figure 2, Multi- Variable Controller Logical Detail 200 shows detail in 
Multi-Variable Controller Logic 1 12. Human Interface Logic 202 provides data output for 
viewing on the screen of Monitor and Keyboard 110 and data input from the keyboard and 
affiliated input devices (for example, trackball, mouse) from Monitor and Keyboard 110. 
Details in Models 204 are shown in Time-Dependent Non-Integrating Characterization 

15 500, Multiple Non- Integrating Model Characterization 600, Non-Integrating Model 

Inertial Characteristics 700, Time-Dependent Integrating Functional Characterization 800, 
Multiple Integrating Model Characterization 900, and Integrating Model Inertial 
Characteristics 1 000 of Figures 5-10. MV Determination Logic 206 executes the 
responsive control logic to determine Manipulated Variable 212 in the context of 

20 Controlled Variable 2 1 4 and Feed Forward Variable 2 1 6 values. Adaptation Logic 208 
provides logic for adapting the control-related values within Models 204 and MV 
Determination Logic 206. Adaptation Logic 208 has Divergent Response Logic 218, Gain 
Determination Logic 220, and Archival Logic 222 logical portions. Real-time Interface 
Logic 210 reads Feed Forward Variable 216 and Manipulated Variable 212 variables from 

25 Control Computer Logic 120 (from data communication enabled by Communication 

Interface 1 06) and communicates these variables to MV Determination Logic 206. Real- 
time Interface Logic 210 also reads Manipulated Variable 212 variables from MV 
Determination Logic 206 and communicates these variables to Control Computer Logic 
120 (through data communication enabled by Communication Interface 106). 

30 

Continuing with attention now to Figure 3, Basic Control System Block Diagram 300 
shows a traditional control system diagram for a controller scheme incorporating a model. 
Controller 302 is the controller (frequently determined as an effective inversion of Model 
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306). Controlled Apparatus 304 is the controller block representation of Operated 
Apparatus 102. Model 306 is a computer-implemented model of aspects of Controlled 
Apparatus 304 under control. Disturbance/CV Summing Point 308 is a summing point for 
Disturbance 312 and output (Controlled Variable) from Controlled Apparatus 304. 
5 Model/CV Summing Point 3 1 0 is a summing point for Model 306 output (a model- 
determined value of a Controlled Variable) and output from Disturbance/CV Summing 
Point 308 (Controlled Variable actuaL value including the influence of any Disturbance 
312). Disturbance 3 12 is any influence on a Controlled Variable not established by 
Controller 402. Setpoint 3 14 represents determinations of either (a) an operating 
10 technician or (b) an automated setpoint determination system as to the desired value of a 
Controlled Variable. 

In Figure 4, Multi- Variable Controller Control System Block Diagram 400 shows a 
control system diagram for a controller scheme incorporating a model scheme according 

15 to the preferred embodiments. Controller 402 is closely affiliated with Controller 302 
respective to abstracted functionality in Multi-Variable Controller Control System Block 
Diagram 400. Controller 402 incorporates Linear Program 426. Controlled Apparatus 
304, Disturbance/CV Summing Point 308, and Disturbance 3 12 are reprised from 
Controlled Apparatus 304. Further detail in Controller 402 is described in Controller 

20 Operation 1700 of Figure 17. Process Disturbance Value Summing Point 404 is a 
summing point for Modeling Error Calculation Block 406 output and Disturbance/CV 
Summing Point 308 output. Output from Process Disturbance Value Summing Point 404 
is input to Controller 402. Modeling Error Calculation Block 406 determines the degree 
of error in model-determination of a Controlled Variable respective to the separate model- 

25 determination instances of the Controlled Variable through use of Fast Model 410, 

Reference Model 408, and Slow Model 412. Reference Model 408 is abstractly equivalent 
to Model 306 except that the model of Reference Model 408 can be modified in use via 
input from Adaptation Block 416. Models of Fast Model 410 and Slow Model 412 are 
determined from the model of Reference Model 408 as further described in Time- 

30 Dependent Non-Integrating Characterization 500, Multiple Non-Integrating Model 
Characterization 600, Non-Integrating Model Inertial Characteristics 700, Time- 
Dependent Integrating Functional Characterization 800, Multiple Integrating Model 
Characterization 900, and Integrating Model Inertial Characteristics 1000 of Figures 5 - 
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10 and in Model Construction Step 1200, Dead-Time Time-Axis Attribute Determination 
1300, Steady-State Time-Axis Attribute Determination 1400, Ramp Portion Determination 
Detail 1500, and Curvalinear Portion Determination Detail 1600 of Figures 12 — 16. 

5 History Block 414 shows Archival Logic 222 in Multi-Variable Controller Control System 
Block Diagram 400. History Block 414 is in data read communication with (a) Controlled 
Variables sent to Controlled Apparatus 304, (b) Manipulated Variables from Controller 
402, and (c) output-model-determined values of the Controlled Variable from Fast Model 
410, Reference Model 408, and Slow Model 412; data received from these control block 

10 elements is archived for use by Adaptation Block 41 6 in adapting Fast Model 410, 

Reference Model 408, Slow Model 412 and Controller 402 as needed. Adaptation Block 
416 is in bilateral data communication with History Block 414 to use the archived values 
in adapting Fast Model 410, Reference Model 408, Slow Model 412, and Controller 402 
as needed. Divergent Response Routine 418 is a logical portion of Adaptation Block 416 

15 for determining control gains for use in Controller 402, Reference Model 408, Fast Model 
410, and Slow Model 412. Gain Determination Routine 420 is a logical portion of 
Adaptation Block 416 for determining if candidate modification data for use in Controller 
402, Reference Model 408, Fast Model 410, and Slow Model 412 would configure the 
system according to Multi- Variable Controller Control System Block Diagram 400 to 

20 provide undesirable divergent control of Operated Apparatus 102. Adaptation Block 416 
inputs adaptive information to Fast Model 410, Reference Model 408, Slow Model 412, 
and Controller 402. In one embodiment, Controller-Adapter Data Path 422 facilitates 
operation of Adaptation Block 416 logic under the direction of Controller 402, and 
Adaptation Block 416 provides inputs via Adapter-Controller Data Path 424. In an 

25 alternative embodiment, Adaptation Block 416 operates in parallel with Controller 402 
without benefit of Controller 402 direction via Controller- Adapter Data Path 422 and 
provides input to Controller 402 via Adapter-Controller Data Path 424. 

Turning now to Figure 5, Time-Dependent Non-Integrating Characterization 500 is useful 
30 in understanding the preferred embodiment modeling approach and shows a depiction of a 
traditional non-integrating response model (a Primary Model) of a Controlled Variable 
change respective to a unit step change in a respective Manipulated Variable such as used 
in defining Model 306. Time-Dependent Non-Integrating Characterization 500 is also the 
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depiction of the response model used in defining Reference Model 408. The Set Of 
Discrete Intervals Of Time 502 are shown on Time Axis 506 with each Discrete Interval 
Of Time 504 being of identical time duration. 

5 An axis for Controlled Variable Magnitude 508 is depicted in orthogonal intersection to 
Time Axis 506 and shows Zero-Time Time- Axis Attribute 512, at the moment (time axis 
zero value) of a unit step change in the respective Manipulated Variable. The function 
between Zero-Time Time-Axis Attribute 512 and Primary Model Dead-Time Time-Axis 
Attribute 516 (Axis 506 locator reference C in Figures 5-10) is usually termed "Dead 

10 Time" in the response model. The model initiates response according to Primary Model 
Curvilinear Portion 522 at Primary Model Dead-Time Time-Axis Attribute 516 (Axis 506 
locator reference C) which continues in time until Primary Model Steady-State Time-Axis 
Attribute 520 (Axis 506 locator reference F in Figures 5-8). From Primary Model Steady- 
State Time-Axis Attribute 520 to Maximum-Time Time-Axis Attribute 514 (Axis 506 

15 locator reference B in Figures 5-10), the Model shows constant response in Primary Model 
Homaloidal Portion 5 1 0 after achieving the modification in the Controlled Variable 
according to Response Gain Attribute 518. 

Turning now to Figure 6, Multiple Non-Integrating Model Characterization 600 reprises 
20 Discrete Intervals Of Time 504, Time Axis 506, Controlled Variable Magnitude 508, 

Primary Model Homaloidal Portion 510, Zero-Time Time-Axis Attribute 512, Maximum- 
Time Time- Axis Attribute 514 (Axis 506 locator reference B), Primary Model Dead-Time 
Time-Axis Attribute 516 (Axis 506 locator reference C), Response Gain Attribute 518, 
Primary Model Steady-State Time-Axis Attribute 520 (Axis 506 locator reference F), and 
25 Primary Model Curvilinear Portion 522 and then shows a Fast (FAST) Model according to 
a function defined by Zero-Time Time-Axis Attribute 512 to Fast Model Dead-Time 
Time- Axis Attribute 602 (Axis 506 locator reference C) to Fast Model Steady-State 
Time- Axis Attribute 604 (Axis 506 locator reference F') to Maximum-Time Time-Axis 
Attribute 514 (Axis 506 locator reference B) having Fast Model Curvilinear Portion 606 
30 (disposed between Fast Model Dead-Time Time-Axis Attribute 602 and Fast Model 

Steady-State Time-Axis Attribute 604) and Fast Model Homaloidal Portion 614 (disposed 
between Fast Model Steady-State Time- Axis Attribute 604 and Maximum-Time Time- 
Axis Attribute 5 14). A Slow SLOW) Model is also shown according to a function defined 
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by Zero-Time Time-Axis Attribute 512 to Slow Model Dead-Time Time-Axis Attribute 
612 (Axis 506 locator reference C") to Slow Model Steady-State Time-Axis Attribute 608 
(Axis 506 locator reference F") to Maximum-Time Time- Axis Attribute 514 having Slow 
Model Curvilinear Portion 610 (disposed between Slow Model Dead-Time Time-Axis 
5 Attribute 612 and Slow Model Steady-State Time-Axis Attribute 608) and Slow Model 
Homaloidal Portion 616 (disposed between Slow Model Steady-State Time-Axis Attribute 
608 and Maximum-Time Time- Axis Attribute 514). The Slow (SLOW) Model and the 
Fast (FAST) Model are both derived from the Primary (REFERENCE) Model according 
to Time-Dependent Non-Integrating Characterization 500. 

10 

Continuing with attention now to Figure 7, Non-Integrating Model Inertial Characteristics 
700 shows Multiple Non-Integrating Model Characterization 600 with further 
identification of Fast Dynamic Response Inertial Characteristic 702 between the Fast 
Model and the Primary Model and of Slow Dynamic Response Inertial Characteristic 704 
15 between the Primary Model and the Slow Model. 

Turning now to Figure 8, a Primary (REFERENCE) Model according to Time-Dependent 
Integrating Functional Characterization 800 shows a depiction of a traditional integrating 
response model (a Primary Model) of a Controlled Variable change respective to a unit 

20 step change in a respective Manipulated Variable such as used in defining Model 306. 
Time Axis 506, Discrete Interval Of Time 504, Zero-Time Time-Axis Attribute 512, 
Maximum-Time Time-Axis Attribute 5 14 (along with respective Time Axis 506 reference 
locator B), and Primary Model Dead-Time Time-Axis Attribute 516 (along with 
respective Time Axis 506 reference locator C) are reprised to show dimensional 

25 alignment with the reference-frame of Time-Dependent Non-Integrating Characterization 
500. Primary Model Ramped Portion 802 differentiates the Primary Model from that 
depicted in Time-Dependent Non-Integrating Characterization 500 and has an attribute of 
a Ramp Rate as the slope of the response as a function of time. 

30 In Figure 9, Multiple Integrating Model Characterization 900 shows the integrating 
response model equivalents of a Fast (FAST) Model in the function defined by points 
Zero-Time Time- Axis Attribute 512, Fast Model Dead-Time Time-Axis Attribute 906 
(Time Axis 506 reference locator C), and Maximum-Time Time- Axis Attribute 514 and a 
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Slow (SLOW) Model in the function defined by Zero-Time Time-Axis Attribute 512, 
Slow Model Dead-Time Time- Axis Attribute 908 (Time Axis 506 reference locator C")> 
and Maximum-Time Time- Axis Attribute 514 (note that Maximum-Time Time- Axis 
Attribute 514 is a singularity "point" respective to point B in all models Time-Dependent 

5 Non-Integrating Characterization 500, Multiple Non-Integrating Model Characterization 
600, Non-Integrating Model Inertial Characteristics 700, Time-Dependent Integrating 
Functional Characterization 800, Multiple Integrating Model Characterization 900, and 
Integrating Model Inertial Characteristics 1000). Primary Model Ramped Portion 802 is 
reprised from the Primary Model along with reference points Discrete Interval Of Time 

10 504, Time Axis 506, Zero-Time Time-Axis Attribute 512, Maximum-Time Time-Axis 
Attribute 514, and Primary Model Dead-Time Time-Axis Attribute 516 (Time Axis 506 
reference locator C). The Fast Model has Fast Model Ramped Portion 902 defined by the 
functional portion connected to Fast Model Dead-Time Time- Axis Attribute 906 and 
Maximum-Time Time-Axis Attribute 5 14. The Slow Model has Slow Model Ramped 

15 Portion 904 defined by the functional portion connected to Slow Model Dead-Time Time- 
Axis Attribute 908 and (at the singularity point denoted by Time Axis 506 reference 
locator B) Maximum-Time Time-Axis Attribute 514. 

Turning now to Figure 10, Integrating Model Inertial Characteristics 1000 shows Multiple 
20 Integrating Model Characterization 900 with further identification of Fast Dynamic 
Response Inertial Characteristic 1002 between the Fast (FAST) Model and the Primary 
(REFERENCE) Model and of Slow Dynamic Response Inertial Characteristic 1004 
between the Primary (REFERENCE) Model and the Slow (SLOW) Model. 

25 Continuing with attention now to Figure 1 1, a consideration of details in methods 

respective to the preferred embodiments, Multi- Variable Controller General Deployment 
Actions 1 1 00 presents a broad overview of steps in using the Model Variable Controller. 
In Model Configuration Step 1 102, Operated Apparatus 102 is operated and data is 
acquired to build the Primary Model and the Primary Model is constructed for use in 

30 Reference Model 408. A DMC (Dynamic Matrix Controller) model is characteristic of the 
initial Primary Model. The Fast and Slow Models of Fast Model 410 and Slow Model 412 
(and Multiple Non-Integrating Model Characterization 600, Non-Integrating Model 
Inertial Characteristics 700, Multiple Integrating Model Characterization 900, and 
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Integrating Model Inertial Characteristics 1000) are also built according to further detail of 
Model Construction Step 1200, Dead-Time Time- Axis Attribute Determination 1300, 
Steady-State Time-Axis Attribute Determination 1400, Ramp Portion Determination 
Detail 1500, and Curvalinear Portion Determination Detail 1600 of Figures 12-16. The 
5 Primary Model is also inverted to provide the response of Controller 402. In Initial 
Tuning Step 1 104, Multi- Variable Controller Control System Block Diagram 400 is 
further configured as a control engineer defines tuning parameters. In Operational Step 
1 106, Multi- Variable Controller Control System Block Diagram 400 is used to control 
Operated Apparatus 102 according to Controller Operation 1700 of Figures 17A and 17B; 

10 in executing this step, disturbance estimation operates within the logic of MV 

Determination Logic 206 to establish, from the set of models, an estimated modeling error 
value in interaction of one Controlled Variable with all the Manipulated Variables. An 
estimated process disturbance value from the Controlled Variable magnitude and the 
estimated modeling error value is also calculated. The estimated model error value and 

15 the estimated process disturbance error value are then used to determine the need for 
further adaptation of either tuning data or model data. 

In providing a discussion of the methods used in the described embodiments, Figures 12 to 
20 present useful sets of key virtual steps showing ordered progressions within the general 

20 motif of a flowchart; as should be apparent however, especially in view of certain parallel 
and alternative paths which are presented without protocol for branching, these flowcharts 
are not meant to rigorously depict a design level template for computer logic coding. A 
conversion from the methods and processes depicted in the figures to a coding flowsheet 
and thence to executable code should be apparent and straightforward given the figures 

25 and discussion. 

Turning now to Figure 12, Model Construction Step 1200 shows a first level of Model 
Configuration Step 1 102 detail in the modification of each Primary Model to a respective 
Fast Model and to a respective Slow Model. Fast Model Curvalinear Portion Modification 
30 Step 1 202 initiates the path for Fast Model Curvilinear Portion 606 definition; and Slow 
Model Curvalinear Portion Modification Step 1204 initiates the path for Slow Model 
Curvilinear Portion 610 definition. 
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Figures 13 — 16 show further details in specific considerations of the models and the 
overview depiction of Model Construction Step 1200. Dead-Time Time-Axis Attribute 
Determination 1300 of Figure 13 shows detail in Primary Model Dead-Time Time- Axis 
Attribute 516 definition. Steady-State Time-Axis Attribute Determination 1400 of Figure 

5 14 shows detail in Primary Model Steady-State Time- Axis Attribute 520 definition. Ramp 
Portion Determination Detail 1500 of Figure 1 5 shows detail in Primary Model Dead- 
Time Time-Axis Attribute 516 and Primary Model Ramped Portion 802 attribute 
determination in Time-Dependent Integrating Functional Characterization 800 and 
Multiple Integrating Model Characterization 900. Curvalinear Portion Determination 

10 Detail 1600 of Figure 16 shows detail in Fast Model Curvilinear Portion 606 and Slow 
Model Curvilinear Portion 610 determination from Primary Model Curvilinear Portion 
522. 

Turning now to Figure 17, Controller Operation 1700 shows the operational process of 

15 MV Determination Logic 206 and Controller 402. In History Update Step 1702, data for 
existing MV, FF, and CV variables is transmitted to History Block 414 for archival and 
use in Adaptation Block 416. The history is built for N input (CV) variables, L feed- 
forward (FF) variables, and M output (MV) variables. In Adaptation Decision Step 1704, 
an estimated process disturbance value (from the Controlled Variable magnitude) and the 

20 estimated modeling error value are calculated. The estimated model error value and the 
estimated process disturbance error value are then used to determine the need for further 
adaptation of either tuning data or model data. Details in Model Adaptation Step 1722 
(given a YES answer from Adaptation Decision Step 1704) are discussed in Adaptation 
Methodology Detail 1800 of Figures 18A- 18E. In CV Prediction Step 1706, the Primary 

25 (Reference), Slow, and Fast Models are used to predict steady-state Controlled Variable 
values. The predictions are done for N input (CV) variables, L feed-forward (FF) 
variables, and M output (MV) variables. Further detail in this is shown in Future CV 
Requirement Definition Detail 1900 of Figure 19. In Steady-State MV Definition Step 
1708, Linear Program 426 is called to define steady-state Manipulated Variable values. 

30 These values are defined for N input (CV) variables and M output (MV) variables. In 

0 

Dynamic Matrix Build Step 1710, the Dynamic (ATA) Matrix is rebuilt if tuning is to be 
changed, if the models are to be changed, or if this is the first execution instance of the 
process of Controller 402. The ATA Matrix has a dimension of M x M, where M = N (CV 
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variables) multiplied by the number of future MV moves for each MV. In Future CV 
Requirement Definition Step 1712, necessary future shifts in Controlled Variable values 
are determined from setpoints and other predicted future values as acquired from the 
database of Control Computer Logic 120 in CV Data Acquisition Step 1720. Further 

5 detail in Future CV Requirement Definition Step 171 2 is shown in CV Prediction Detail 
2000 of Figure 20. The shifts are determined for M output variables. In MV Change 
Definition Step 1714, the ATE matrix and Dynamic Matrix are solved to define 
incremental changes in Manipulated Variables. In MV Implementation Step 1716, the 
incremental Manipulated Variable changes are implemented and the process returns to 

10 Data Acquisition Step 1718; this affects each of N input variables. In Data Acquisition 
Step 1718, MV, CV, and FF variables are read from Control Computer Logic 120 (via 
Communication Interface 106). In CV Data Acquisition Step 1720, setpoints and other 
predicted future values are acquired from the database of Control Computer Logic 120 
along with necessary future shifts in Controlled Variable values. 

15 

In Figure 18, Adaptation Methodology Detail 1800 shows specifics in Model Adaptation 
Step 1722 (given a YES answer from Adaptation Decision Step 1704). In Forgetting 
Factor Modification Decision 1804, the "forgetting factor" is applied to the most senescent 
data in Archival Logic 222 via use in determining individual discount factors. In C Matrix 
20 Build Step 1 802 and other steps of Adaptation Methodology Detail 1 800, "differential 
form" is sometimes also termed "difference form". The C Matrix is dimensioned as M x 
N where N is the number of MV and FF variables and M is the number of increments in 
time horizon used. 

25 Continuing with attention now to Figure 1 9, Future CV Requirement Definition Detail 
1900 shows process detail in Future CV Requirement Definition Step 1712 of Controller 
Operation 1700. In Alternative CV Error Determination Step 1902, an alternative error 
value is determined for the current and predicted CV values. In example: 

30 Equations 58 - 60 

L = maximum predicted CV - minimum predicted CV 
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X = current CV/L 
error = x * (1.0-exp(-x*x)) 
5 OR 

eiror = x * (1.0-L0/(l+a*x A n)) 

where "a" is defined to enable robust tuning (note that "a" is, in one embodiment, a 
10 function of constraint proximity). 

In Figure 20, CV Prediction Detail 2000 presents further detail in Future CV Requirement 
Definition Step 1712. In CV Error Correction Step 2002, each error for each Controlled 
Variable is evaluated via the relationship: 

15 

Equation 61 

Error = setpoint — function(Fast predicted CV, Slow predicted CV) 

20 where the function is a straightforward arithmetic mean or, alternatively, the function is 
based upon constraint proximity in the CV. 

The benefit of the multi- variable controller when compared to a traditional DMC 
controller is demonstrated in the following Examples from simulation in conjunction with 
25 Figures 21 to 26. 

In the example, a simulation of an operating apparatus is controlled in Figures 21, 23, and 
25 with a classic DMC. In respectively comparative Figures 22, 24, and 26, the same 
simulation of the operating apparatus is controlled using the robust multi- variable 
30 controller of the preferred embodiments. In the Figure 21-22 comparison, the model is 
"affected" with an identical shift in a model parameter between the simulation of the plant 
and the model used in the controller. In the Figure 23-24 comparison, the model is 
"affected" with an identical shift in a gain parameter between the simulation of the plant 

-50- 



WO 02/03150 



PCT/US01/18111 



and the model used in the controller. In the Figures 25-26 comparison, the model is 
"affected" with an identical controller model (characteristic value) inversion. In reading 
the simulation output time charts, note that the right hand scale of the time charts defines 
the quantitative value of the MV variable, whereas the left hand scale defines the SP and 
5 CV quantitative values. 

Figure 21 presents output 2100 from a simulator for a regular DMC operating in a 
situation of model mismatch in modeling parameters. Output 2100 contains time chart 
2102 for the actions of CVi and MVi, time chart 2104 for the actions of CV 2 and MV 2 , 

10 and time chart 21 06 for the actions of CV 3 and MV 3 . Figure 2 1 shows the regular DMC in 
the situation of model mismatch between the simulation of the operating apparatus and the 
model used in the controller (the case is a 3x3 with model mismatch in only the first 
model; in this regard, the dead time is incorrectly modeled). The move suppression factor 
is set to unity. Equal concern errors are also all unity. The simulation is normalized 

15 internally so that all values of CV and MV begin at 50.0. The process gains are usually in 
the 0 to 3.0 range. The figure shows the instability induced in the DMC by a single small 
error in the dead time model in only one model out of 9 (that is, 3x3). Each other model 
parameter is strictly equal between the apparatus simulation and the model of the 
controller. The dead time error in the first model is about 20 percent. This error makes the 

20 controller unstable in all the 3 CV's; the traditional cure is to increase move suppression 
but at the expense of controller reaction time to external disturbance and to set point 
changes. 

Figure 22 presents output 2200 from a simulator for the robust controller of the preferred 
25 embodiments operating in the situation of model mismatch of Figure 21 . Output 2200 
contains time chart 2202 for the actions of CVj and MVi, time chart 2204 for the actions 
of CV 2 and MV 2 , and time chart 2206 for the actions of CV 3 and MV 3 . Figure 22 shows 
the robust controller in face of the same model mismatch as Figure 21 (the case is the 
same 3x3 with model mismatch in only the first model; the dead time is again incorrectly 
30 modelized). The move suppression factor is set to unity. Equal concern errors are also all 
unity. The simulation is normalized internally so that all values of CV and MV start at 
50.0. The figure shows the resulting stability of the robust multivariable controller. The 
controller is not gaining stability at the expense of the reaction time, and, therefore, 
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disturbances and set point changes can be handled faster than the regular DMC. The 
robustness is so substantial that move suppression can be zero (i.e turned off) if there is 
any reason to do so. 

5 Figure 23 presents output 2300 from a simulator for a regular DMC operating in a 
situation of model mismatch in gains. Output 2300 contains time chart 2302 for the 
actions of CVi and MVi, time chart 2304 for the actions of CV 2 and MV 2 , and time chart 
2306 for the actions of CV 3 and MV3. Figure 23 shows the regular DMC without 
adaptation. The case is a 3x3 with model mismatch in only the gains of all models. The 
10 move suppression factor is set to unity. Equal concern errors are also all unity. The 
simulation is normalized internally so that all values of CV and MV start at 50.0. The 
process gains are usually in the 0 to 3.0 range with errors in the range 0 to 50 percent. The 
figure shows the instability induced in the DMC by the gain errors. The DMC controller 
is unstable in all the 3 CV's. 

15 

Figure 24 presents output 2400 from a simulator for the robust controller of the preferred 
embodiments operating in the situation of model mismatch of Figure 23. Output 2400 
contains time chart 2402 for the actions of CVi and MVi , time chart 2404 for the actions 
of CV 2 and MV 2 , and time chart 2406 for the actions of CV3 and MV 3 . Figure 24 shows 

20 the robust controller with adaptation. The case is the same 3x3 as used in Figure 23 with 
model mismatch in only the gains of all models. The move suppression factor is set to 
unity. Equal concern errors are also all unity. The simulation is normalized internally so 
that all values of CV and MV start at 50.0. The process gains are usually in the 0 to 3.0 
range with errors in the range 0 to 50 percent. The figure shows some initial instability in 

25 the robust controller (area of 2408); this is a learning period. After this learning period, 
the controller exhibits nearly perfect response to set point changes since it derives, from 
past data, the correct model gains that match the actual process to the control model. Note 
also, in comparing the peak in the area of 2408 of Figure 24 with the comparable peak 
area of 2308 of Figure 23, that the controller of Figure 24 shows less overshoot above the 

30 setpoint SP; this demonstrates the efficiency with which the described embodiment of the 
multi-model controller reacts to the operating system. 



-52- 



WO 02/03150 



PCT/US01/18111 



Figure 25 presents output 2500 from a simulator for a regular DMC operating in a 
situation of controller model (characteristic value) inversion. Output 2500 contains time 
chart 2502 for the actions of CVj and MVj, time chart 2504 for the actions of CV 2 and 
MV 2 , and time chart 2506 for the actions of CV 3 and MV 3 . Note that (a) the set point 
5 change at time t=l 00 triggers the inverted response and (b) the controller diverges very 
rapidly since the move suppression is set to unity. 

Figure 26 presents output 2600 from a simulator for the robust controller of the preferred 
embodiments operating in the situation of model (characteristic value) inversion of Figure 

10 25. Output 2600 contains time chart 2602 for the actions of CVi and MVi, time chart 2604 
for the actions of CV 2 and MV 2 , and time chart 2606 for the actions of CV 3 and MV 3 . 
Again, the set point change at time t=100 triggers the inverted response, and the controller 
initially diverges - but not as rapidly as the regular DMC of Figure 25. Inversion 
detection then triggers inversion of the characteristic value in all models; this process 

15 effectively occurs at time t= 120. As should be apparent, the inversion is initiated at an 
earlier relative moment when a different choice of controller parameters are used. The 
move suppression factor is set to unity. A comparison of Figures 25 and 26 shows the 
value in the described embodiments of (a) determining divergent response behavior in the 
controlling process and then (b) counteracting the identified divergent response behavior 

20 to stabilize the controller. 

The described embodiments are achievable within a number of computer system 
architectural alternatives. In one alternative, an embodiment is facilitated within the 
context of a multi-process environment wherein different databases, data sections, and 

25 logical engines (logical sub-sections which read data, write data, calculate data, and make 
decisions in data computational processes) within the computer-implemented logic are 
simultaneously installed and activated with dynamically active data transfer linkages, 
facilitated either directly or indirectly via the use of a data common and/or an application 
program interface (APIs). In another alternative, the different databases, data sections, 

30 and logical engines are facilitated within the context of a single process environment 
wherein different components are sequentially activated by an operating technician with 
linkages facilitated either directly or indirectly via the use of data commons or data 
schema dedicated to interim storage. In yet another alternative, the different databases, 
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data sections, and logical engines are deployed within the context of a single process 
environment wherein (a) some components of the different databases, data sections, and 
logical engines are accessed and activated by an operating technician with linkages 
facilitated either directly or indirectly via the use of data commons or data schema 

5 dedicated to interim storage, and (b) the other components within the different databases, 
data sections, and logical engines are accessed and activated by calls from previously 
activated with linkages facilitated either directly or indirectly via the use of data commons 
or data schema dedicated to interim storage. In one alternative, the multi-variable 
controller is implemented and executed on one physical computer. In another alternative, 

10 the controller is facilitated on different platforms where the results generated by one 
engine are transferred by an operating technician to a second or other plurality of the 
different databases, data sections, and logical engines executing on different computer 
platforms, although a separate operating system is needed on each platform. In yet 
another alternative, controller is facilitated on a plurality of computer platforms 

15 interconnected by a computer network, although a separate operating system is needed on 
each platform and the operating system further incorporates any networking logic that is 
needed to facilitate necessary communications via such a computer implemented 
communication network. A relatively small controller according to the described 
embodiments is deployed on a computer having an Intel 80486 CPU with a 33 MHz clock, 

20 10 Megabytes of RAM Memory, and a 100 Megabyte Hard Disk using a Windows '95 
operating system from Microsoft Corporation. A larger controller according to the 
described embodiments is deployed on a Vaxstation 4000m90 having 128 Megabytes of 
RAM and a (at least) 500 Megabyte Hard Disk from Compaq Computer Corporation. 
Many of the different gradations of architectural deployment within the context of the 

25 above overview are considered by the applicants to be generally apparent, and the 

illustration of present invention can be conveniently modified by those of skill, given the 
benefit of this disclosure, to achieve the utility of the present invention within the context 
of the above computer system architectural alternatives without departing from the spirit 
of the present invention once given the benefit of the disclosure. 
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I claim: 

1 . A computer-implemented system for controlling the operation of a processing 
apparatus having at least one independently controlled manipulated variable and at least 

5 one controlled variable responsive to said manipulated variable, characterized by: 

a set of at least two models for defining an expected variation in magnitude for each 
controlled variable as a respective function of each manipulated variable, said model set 
having a dynamic response inertial characteristic; and 

10 

means for implementing a change defined from said model set to modify said manipulated 
variable in said processing apparatus. 

2. The system of Claim 1 further characterized by means for adapting said models 
1 5 during real-time use. 

3. The system of Claim 1 further characterized by means for determining divergent 
response behavior in said controlling. 

20 4. A computer-implemented system for controlling the operation of a processing . 
apparatus having at least one independently controlled manipulated variable and at least 
one controlled variable responsive to at least one said manipulated variable and further 
responsive to process disturbances originating independently of said manipulated variable, 
characterized by: 

25 

means for measuring the magnitude of each controlled variable; 

a set of at least two models for defining an expected variation in magnitude for each 
controlled variable as a respective function of a manipulated variable disturbance instance 
30 in each manipulated variable; 



-55- 



WO 02/03150 



PCT/US01/18111 



means for determining, from said set of models, an estimated modeling error value in 
interaction of one said controlled variable with all said manipulated variables in achieving 
a manipulated variable modification; 

5 means for determining an estimated process disturbance value from said controlled 
variable magnitude and said estimated modeling error value; 

means for defining a first portion of desired change in the present value of at least one 
manipulated variable from said estimated modeling error value; 

10 

means for defining a second portion of desired change in the present value of at least one 
manipulated variable from said estimated process disturbance value; and 

means for implementing said first and second desired change portions to modify said. 
1 5 manipulated variable. 

5. The system of Claim 4 wherein said model set incorporates a dynamic response 
inertial characteristic, said system further characterized by means for acquiring an inertial 
characteristic value so that said dynamic response inertial characteristic can be achieved in 

20 said model set. 

6. The system of Claim 5 further characterized by means for adapting said model set 
during real-time use. 

25 7. The system of Claim 4 further characterized by means for determining divergent 
response behavior in said controlling. 

8. A computer-implemented system for controlling the operation of an apparatus 
having at least one independently controlled manipulated variable and at least one 
30 controlled variable responsive to at least one said independently controlled manipulated 
variable, characterized by: 
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means for defining a set of consecutive discrete intervals of time in a time-dependent 
function; 

means for introducing, in each independently controlled manipulated variable, a 
5 manipulated variable disturbance instance of predefined magnitude, said disturbance 
instance prompting a response in each controlled variable; 

means for measuring the magnitude of each controlled variable; 

10 a controller; 

means for defining, respective to each response, at least one time-dependent functional 
characterization of said controlled variable magnitude over said set of consecutive discrete 
intervals of time on a time-axis, said functional characterization having a zero-time time- 

15 axis attribute, a maximum-time time-axis attribute, a dead-time time-axis attribute, a 
response gain attribute, a ramp-rate attribute, a steady-state time-axis attribute, a 
curvilinear portion disposed between the time-axis position of said dead-time time-axis 
attribute and the time-axis position of said steady-state time-axis attribute, a ramped 
portion disposed between the time-axis position of said dead-time time-axis attribute and 

20 the time-axis position of said maximum-time time-axis attribute, and a homaloidal portion 
disposed between the time-axis position of said steady-state attribute and the time-axis 
position of said maximum-time time-axis attribute, each discrete time interval for one said 
response having the same time duration, said homaloidal portion having a zero value for 
an integrating controlled variable response, said ramped portion having a zero value for a 

25 non-integrating controlled variable response, and each functional characterization for one 
said response having identically valued zero-time time-axis attributes, response gain 
attributes, ramp-rate attributes, and maximum-time time-axis attributes; 

means for acquiring a first said time-dependent functional characterization respective to 
30 the response from measuring an effected change in said magnitude of a controlled variable 
after introducing said disturbance instance, said first functional characterization having a 
first said dead-time time-axis attribute, a first said steady-state time-axis attribute, a first 
said curvilinear portion, a first said homaloidal portion, and a first said ramped portion 
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having its functional derivative equivalent to said ramp-rate attribute at said maximum- 
time time-axis attribute; 

means for inverting said first time-dependent functional characterization into said 
5 controller; 

means for deriving a second said time-dependent functional characterization from said 
first time-dependent functional characterization, said second functional characterization 
having a second dead-time attribute in first predefined diminishing offset from said first 
10 dead-time attribute, a second steady-state attribute in second predefined diminishing offset 
from said first steady-state attribute, a second curvilinear portion in first predefined 
functional offset from said first curvilinear portion, a second homaloidal portion in 
extrapolation of said first homaloidal portion, and a second ramped portion in 
extrapolation of said first ramped portion; 

15 

means for deriving a third said time-dependent functional characterization from said first 
time-dependent functional characterization, said third functional characterization having a 
third dead-time attribute in first predefined superadditive offset from said first dead-time 
attribute, a third steady-state attribute in second predefined superadditive offset from said 
20 first steady-state attribute, a third curvilinear portion in second predefined functional offset 
from said first curvilinear portion, and a third homaloidal portion in truncation of said first 
homaloidal portion, and a third ramped portion in truncation of said first ramped portion; 

means for determining a desired change in the value of a manipulated variable in real-time 
25 from said second time-dependent functional characterization, said third time-dependent 
functional characterization, the magnitude of at least one controlled variable, and said 
controller; and 

means for implementing said desired change to modify said manipulated variable. 

30 

9. The computer-implemented system of Claim 8, said means for determining a 
desired change in the value of a manipulated variable in real-time determining said desired 
change from said first time-dependent functional characterization, said second time- 
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dependent functional characterization, said third time-dependent functional 
characterization, the magnitude of at least one controlled variable, and said controller. 

10. The computer-implemented system of Claim 9 wherein said first, second, and third 
5 time-dependent functional characterizations define three models in a model set, said model 

set incorporating a first dynamic response inertial characteristic between said first and 
second time-dependent functional characterizations and a second dynamic response 
inertial characteristic between said first and third time-dependent functional 
characterizations, said system further characterized by means for acquiring a first inertial 
10 characteristic value and a second inertial characteristic value so that said first and second 
dynamic response inertial characteristics can be respectively achieved in said model set. 

1 1 . The system of Claim 1 0 further characterized by means for adapting said model set 
during real-time use. 

15 

12. A computer-implemented system for controlling the operation of a processing 
apparatus having at least one independently controlled manipulated variable and at least 
one controlled variable responsive to at least one said manipulated variable and further 
responsive to process disturbances originating independently of said manipulated variable, 

20 characterized by: 

means for measuring the magnitude of each controlled variable; 

a model for defining an expected variation in magnitude for each controlled variable as a 
25 respective function of a manipulated variable disturbance instance in each manipulated 
variable; 

means for determining, from said model, an error value in interaction of one said 
controlled variable with all said manipulated variables in achieving a manipulated variable 
30 modification; 

means for defining a desired change in the present value of at least one manipulated 
variable from said error value; 
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means for implementing said desired change to modify said manipulated variable; and 
means for determining divergent response behavior in said controlling. 

5 

1 3 . The computer-implemented system of Claim 1 2 further characterized by means for 
counteracting said divergent response behavior. 

14. A method for computer-implemented controlling of a processing apparatus having 
10 at least one independently controlled manipulated variable and at least one controlled 

variable responsive to said manipulated variable, characterized by the steps of: 

defining an expected variation in magnitude for each controlled variable as a respective 
function of each manipulated variable from a set of at least two models, said model set 
15 having a dynamic response inertial characteristic; and 

implementing a change defined from said model set to modify said manipulated variable 
in said processing apparatus. 

20 15. The method of Claim 14 further characterized by the step of adapting said models 
during real-time use. 

16. The method of Claim 14 further characterized by the step of determining divergent 
response behavior in said controlling. 

25 

17. A method for computer-implemented controlling of a processing apparatus having 
at least one independently controlled manipulated variable and at least one controlled 
variable responsive to at least one said manipulated variable and further responsive to 
process disturbances originating independently of said manipulated variable, characterized 

30 by the steps of: 

measuring the magnitude of each controlled variable; 
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defining an expected variation in magnitude for each controlled variable as a respective 
function of a manipulated variable disturbance instance in each manipulated variable from 
a set of at least two models; 

5 determining, from said set of models, an estimated modeling error value in interaction of 
one said controlled variable with all said manipulated variables in achieving a manipulated 
variable modification; 

determining an estimated process disturbance value from said controlled variable 
10 magnitude and said estimated modeling error value; 

defining a first portion of desired change in the present value of at least one manipulated 
variable from said estimated modeling error value; 

15 defining a second portion of desired change in the present value of at least one 
manipulated variable from said estimated process disturbance value; and 

implementing said first and second desired change portions to modify said manipulated 
variable. 

20 

18. The method of Claim 17 wherein said model set incorporates a dynamic response 
inertial characteristic, said method further characterized by the step of acquiring an inertial 
characteristic value so that said dynamic response inertial characteristic can be achieved in 
said model set. 

25 

1 9. The method of Claim 1 8 further characterized by the step of adapting said model 
set during real-time use. 

20. The method of Claim 17 further characterized by the step of determining divergent 
30 response behavior in said controlling. 



21 . A method for computer-implemented controlling of an apparatus having at least 
one independently controlled manipulated variable and at least one controlled variable 
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responsive to at least one said independently controlled manipulated variable, 
characterized by the steps of: 

defining a set of consecutive discrete intervals of time in a time-dependent function; 

5 

introducing, in each independently controlled manipulated variable, a manipulated 
variable disturbance instance of predefined magnitude, said disturbance instance 
prompting a response in each controlled variable; 

10 measuring the magnitude of each controlled variable; 

providing a controller; 

defining, respective to each response, at least one time-dependent functional 
15 characterization of said controlled variable magnitude over said set of consecutive discrete 
intervals of time on a time-axis, said functional characterization having a zero-time time- 
axis attribute, a maximum-time time-axis attribute, a dead-time time-axis attribute, a 
response gain attribute, a ramp-rate attribute, a steady-state time-axis attribute, a 
curvilinear portion disposed between the time-axis position of said dead-time time-axis 
20 attribute and the time-axis position of said steady-state time-axis attribute, a ramped 

portion disposed between the time-axis position of said dead-time time-axis attribute and 
the time-axis position of said maximum-time time-axis attribute, and a homaloidal portion 
disposed between the time-axis position of said steady-state attribute and the time-axis 
position of said maximum-time time-axis attribute, each discrete time interval for one said 
25 response having the same time duration, said homaloidal portion having a zero value for 
an integrating controlled variable response, said ramped portion having a zero value for a 
non-integrating controlled variable response, and each functional characterization for one 
said response having identically valued zero-time time-axis attributes, response gain 
attributes, ramp-rate attributes, and maximum-time time-axis attributes; 

30 

acquiring a first said time-dependent functional characterization respective to the response 
from measuring an effected change in said magnitude of a controlled variable after 
introducing said disturbance instance, said first functional characterization having a first 
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said dead-time time-axis attribute, a first said steady-state time-axis attribute, a first said 
curvilinear portion, a first said homaloidal portion, and a first said ramped portion having 
its functional derivative equivalent to said ramp-rate attribute at said maximum-time time- 
axis attribute; 

5 

inverting said first time-dependent functional characterization into said controller; 

deriving a second said time-dependent functional characterization from said first time- 
dependent functional characterization, said second functional characterization having a 

10 second dead-time attribute in first predefined diminishing offset from said first dead-time 
attribute, a second steady-state attribute in second predefined diminishing offset from said 
first steady-state attribute, a second curvilinear portion in first predefined functional offset 
from said first curvilinear portion, a second homaloidal portion in extrapolation of said 
first homaloidal portion, and a second ramped portion in extrapolation of said first ramped 

15 portion; 

deriving a third said time-dependent functional characterization from said first time- 
dependent functional characterization, said third functional characterization having, a third 
dead-time attribute in first predefined superadditive offset from said first dead-time: 
20 attribute, a third steady-state attribute in second predefined superadditive offset from said 
first steady-state attribute, a third curvilinear portion in second predefined functional offset 
from said first curvilinear portion, and a third homaloidal portion in truncation of said first 
homaloidal portion, and a third ramped portion in truncation of said first ramped portion; 

25 determining a desired change in the value of a manipulated variable in real-time from said 
second time-dependent functional characterization, said third time-dependent functional 
characterization, the magnitude of at least one controlled variable, and said controller; and 

implementing said desired change to modify said manipulated variable. 

30 

22. The method of Claim 21, said step of determining a desired change in the value of a 
manipulated variable in real-time determining said desired change from said first time- 
dependent functional characterization, said second time-dependent functional 
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characterization, said third time-dependent functional characterization, the magnitude of at 
least one controlled variable, and said controller. 

23. The method of Claim 22 wherein said first, second, and third time-dependent 
5 functional characterizations define three models in a model set, said model set 

incorporating a first dynamic response inertial characteristic between said first and second 
time-dependent functional characterizations and a second dynamic response inertial 
characteristic between said first and third time-dependent functional characterizations, said 
method further characterized by the step of acquiring a first inertial characteristic value 
10 and a second inertial characteristic value so that said first and second dynamic response 
inertial characteristics can be respectively achieved in said model set. 

24. The method of Claim 23 further characterized by the step of adapting said model 
set during real-time use. 

15 

25. A method for computer-implemented controlling of a processing apparatus having 
at least one independently controlled manipulated variable and at least one controlled 
variable responsive to at least one said manipulated variable and further responsive to 
process disturbances originating independently of said manipulated variable, characterized 

20 by the steps of: 

measuring the magnitude of each controlled variable; 

a model for defining an expected variation in magnitude for each controlled variable as a 
25 respective function of a manipulated variable disturbance instance in each manipulated 
variable; 

determining, from said model, an error value in interaction of one said controlled variable 
with all said manipulated variables in achieving a manipulated variable modification; 

30 

defining a desired change in the present value of at least one manipulated variable from 
said error value; 
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implementing said desired change to modify said manipulated variable; and 

determining divergent response behavior in said controlling. 

5 26. The method of Claim 25 further characterized by the step of counteracting said 
divergent response behavior. 
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FIG. 14 

I400 



C^_startJ^> 



Starting with last two values of 
model 




YES 



Steady State Time is at N+ 1 
Steady State Value = i (at time=N+l) 
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FIG. 15 

I500 

C^START^ 



Read last two values of model 
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FIG. 16 



I600 



Curvalinear portion estimation 




Starts at DEAD TIME and 
Ends at STEADY STATE 




Select indice starting with first 
curvalinear value in the reference model 



Derive j as an expanded indice 
j = i*N/M 




Fast model(i)= model(j+l)*FC + 
model(j)*(10-FC) 




■r 



Nexti 




(Calculate new span 
from Robust factor 



/ Define M such that 
I M>N 



Select indice starting with first curvalinear 
value in the reference model 



Derive j as a contracted indice 
j = i*N/M 




Slowmodel(i)=model(j+l)*FC +\ 
model(j)*(1.0-FC) J 



Nexti 
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1700 
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FIG. 17B 



o 



1700 

/ 



Use model and MVs and FFs to predict steady state CVs: 
-with FAST model 

-with REFERENCE model (optional) 
-with SLOW model 



1706 



Call Linear Program to establish steady-state MVs and (optionally) 
predicted steady-state CVs 



-1708 



Build Dynamic (ATA) Matrix 
(if tuning changed or 
models changed or 

if this is the first execution instance of this method) 



--I7IO 



Define required future changes in 
CVs from shifted predictions of 
controlled CVs 



1712 



Read setpoints and their 
predicted future values 
respective to the CVs from 
the Control Computer 



1720 



Define ATE matrix and solve Dynamic Matrix to 
derive MV incremental changes 



1714 



Implement MV incremental 
chang s: 

New M Vs= old MVs + MV 
incremental change 



1716 
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FIG. 18B 



I800 




Build C matrix: 

each element is the sum of the model elements (in differential form) 

multiplied by 
changes in MV (i.e. MV in differential form) 



-I802 



Expand C matrix for FF variables, 
each element is the sum of the model elements (in differential form) 

multiplied by 
changes in FF (i.e. FF in differential form) 



Apply discount factors to C matrix. 
Factors are applied on each column 
(inner product by row) 



Build CtC matrix 



Build CtP matrix 



I 



Solve CtC = CtP 
to yield : A Gains factor matrix 
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1 800 




From REFERENCE model: 
Evaluate each characteristic value 
(eigenvalue) model M and retain the vector 



For new model: 

Evaluate each characteristic value 
(eigenvalue) model M' 




For each characteristic value 
(eigenvalue) in each model that is 
greater than e (the threshold) 
retain the new space 
transformation vector 
corresponding to the accepted 
characteristic value (eigenvalue) 
to update the gain matrix 




Retain newly calculated characteristic values 
(eigenvalues) in history area 



i 



For each characteristic value (eigenvalue), 
filter and estimate current value using data 
from the historical archive. 
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Express last change in each MV as 
a transformed MV according to T 





f ' ^ 
Express last change in each CV as a 

transformed C V according to the 

inverse of T 




1 


r 



For each transformed MV, derive the 
corresponding characteristic value 
(eigenvalue) by dividing any change in MV (if 
non-zero and greater than a validation number) 



I 





Retain history of the newly 
calculated characteristic 
values (eigenvalues) 




1 


f 



>v For each characteristic value - 

(eigenvalue), filter and estimate 
current value using historical data 



M800 
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Initialize 


/ 


1 


r 





For current time in prediction, find the min and max of each predicted CV value, 
(for both SLOW and FAST prediction) 

Optionally include other predicted CV values (example: REFERENCE model) 



Set error to zero 




NO 



Error = current CV value - 

maximum predicted CV value 



1902 



(alternative) 

Error = function( current CV, 
minimum predicted CV, maximum 
predicted CV) 
V J 



Error = current C V value - 

minimum predicted CV 



Increase or decrease all predictions (SLOW, FAST and 
optionally REFERENCE models and any other) by amount 
of calculated error 



1 


r 




End 





WO 02/03150 



25/3 7 



PCT/USO 1/18111 



FIG. 20 



Initialize 




f 


Read controlled variables setpoints and all time limits, 

V J 




r 


Initialize counters for time increments and position in 
vector 


\ 

r 


f 

> 



2000 



Calculate actual required movements in each CV (i.e. the 
set of "errors" to be corrected) 



2002 







Apply a scaling sensitivity factor to each resulting error 




V 



Build ATE from derived vector set 
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